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q_) ' Abstract 

We introduce and carefully study a natural probability measure over the numerical range of a 
complex matrix A E M„(C). This numerical measure [ia can be defined as the law of the random 
variable (AX, X) € C when the vector X 6 C" is uniformly distributed on the unit sphere. If the 
matrix A is normal, we show that \ia has a piecewise polynomial density Ja, which can be identified 
fy | with a multivariate £?-spline. In the general (nonnormal) case, we relate the Radon transform of 

\ia to the spectrum of a family of Hermitian matrices, and we deduce an explicit representation 
formula for the numerical density which is appropriate for theoretical and computational purposes. 
As an application, we show that the density Ja is polynomial in some regions of the complex plane 
which can be characterized geometrically, and we recover some known results about lacunas of 
symmetric hyperbolic systems in 2 + 1 dimensions. Finally, we prove under general assumptions 
that the numerical measure of a matrix A G M n (C) concentrates to a Dirac mass as the size n goes 
■ to infinity. 

(N 
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^ ■ 1 Introduction 

If A € M n (C) is a complex square matrix of size n 6 N*, the numerical range of A is the compact 
^ ■ subset of the complex plane defined by 



X 



W(A) = \{Ax,x) G C 



x £ C n , llxll = 1 



where (x,y) = y*x is the usual scalar product in C n and ||x|| = (x,x) l l 2 . It is quite obvious that 
iy(j4) D c( j 4)i where cr(A) (the spectrum of ^4) is the collection of all eigenvalues of A, and that 
W(j4) = W(U*AU) for any unitary matrix U E U n (C). Moreover, a celebrated result due to Toeplitz 
[21] and Hausdorff [13] asserts that VF(^4) is always a convex subset of the complex plane. In particular, 
iy(j4) contains the convex hull of cr(A), and it is easy to verify that = conv(cr(j4)) if the matrix 

A is normal, namely AA* = A* A. The interested reader is referred to Chapter 1 of |15] for a detailed 
discussion of the various properties of the numerical range, including complete proofs. 

Let dM n = {x 6 C n | ||x|| = 1} be the unit sphere in C n , considered as a real manifold of dimension 
2n — 1. By definition, the numerical range 1T(>1) is the image of the numerical map : dM n — > C 
defined by 

$ A (x) = (Ax,x) , x£dM n . 
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The algebraic and geometric properties of the map $>a have been extensively studied, see [HI [261 HI 
El [161 HZ] - I n particular, the set of all critical values of &a, which we denote by XU C C, has received 
a lot of attention, because this is an interesting object which contains a lot of information on the 
matrix A. For instance, it is known that dW(A) C Ha and W(-A) = com^E^). In addition, there 
exists a real algebraic curve Ca C C ~ M? with the property that E^ = Ca U C' a , where C' A denotes 
the set of all line segments joining pairs of points of Ca at which Ca has the same tangent line |17| . 
Under generic assumptions on A, the bitangent set C' A is empty, and the critical set E^ is therefore 
the union of a finite number of closed curves, one of which is the boundary of the numerical range 
W(A). This distinguished curve is smooth, and encloses all the other ones in its interior. We refer 
to Section [5] below for more details on the geometry of the singular set, and to Section [7] for a few 
concrete examples. 

Our purpose in this paper is to introduce another mathematical quantity which is naturally related 
to the numerical map &a- Given A G M n (C), the numerical measure of A is the probability measure 
HA on C defined by the formula 

[ <f>{z)dn A (z) = I <j>({Ax,x))da(x) , (1) 
./;• J 8M n 

for all continuous functions <j) : C — > C. Here a denotes the Euclidian measure on the unit sphere dM n , 
normalized as a probability measure. In words, the numerical measure is thus the image under the 
numerical map of the normalized Euclidean measure on the unit sphere. Equivalently, if X is a random 
variable that is uniformly distributed on dM n , the numerical measure ha is just the distribution of 
the random variable (AX, X) E C. This probabilistic intepretation will be useful later, especially in 
Section 

Our first goal is to establish a few general properties of the numerical measure ha- It is clear 
by construction that ha is invariant under unitary conjugations of A, namely hu*au = HA for all 
U E U n (C). This is precisely the reason why we used the Euclidean measure on dM n in the definition 
([1]). It is also easy to verify that the support of HA is exactly the numerical range W(-A), see Section [2] 
below. Less obvious, perhaps, is the fact that ha is absolutely continuous with respect to the Lebesgue 
measure A on W(-A), so that we can define the numerical density /a as the Radon-Nikodym derivative 
of ha with respect to A (in the particular situation where VF(^4) reduces to a line segment T, we 
understand A as the one-dimensional Lebesgue measure on T, see Section [2j) We also prove that the 
numerical density /a is strictly positive in the interior of A, a property that can be interpreted as a 
strong version of Hausdorff's theorem [13]. Finally, we shall see that the singular support of ha is 
contained in the critical set T*a, which means that the numerical density /a is smooth outside Y*a- In 
fact, we conjecture that singsupp(/XA) = E^ for all A E M„(C), but this has not been proved yet. 

After these general properties have been established, our next goal is to give a more precise 
description of the numerical density /a- For this purpose, it is convenient to distinguish between 
various cases: 

1. (The scalar case) If W(.A) is reduced to a single point {z}, then A = zl n (where /„ E M n (C) 
denotes the identity matrix) and ha = S z . In this trivial situation, there is of course no need to 
introduce a numerical density. 

2. (The Hermitian case) Assume that n > 2 and that W(A) C C is a line segment. Then we can find 
z G C, 9 G [0, 2tt], and a Hermitian matrix H such that A = zl n + e l6 H. If Ai < A2 < • • • < A n are the 
eigenvalues of H, we shall see in Section [3] that the numerical measure hh is absolutely continuous with 
respect to Lebesgue's measure on R, and that the corresponding density fjj is exactly the normalized 
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B -spline of degree n — 2 with knots Ai, . . . , A n [7]. In particular, fn is polynomial of degree n — 2 on 
each interval [Ai, Aj+i], vanishes identically outside [Ai, A n ], and is continuous at each point Aj together 
with its derivatives up to order di = n — 2 — mi, where m, > 1 is the multiplicity of Aj as an eigenvalue 
of H (if di < 0, then fjj is discontinuous at Aj.) This gives an explicit representation of the numerical 
measure and the measure \ia is the image of [iu under the affine isometry n4z-(- e w. 

3. (The normal case) Suppose now that n > 3 and that A G M n (C) is a normal matrix whose 
spectrum cr(A) is not contained in a line segment. Then W(A) = conv(<r(^4)) is a convex polygon 
with nonempty interior, and it turns out that the numerical density Ja is the bivariate B-spline of 
degree n — 3 whose knots are the eigenvalues of A. Here we refer to the work of W. Dahmen [6] for 
the definition and the main properties of multivariate l?-splines. In this particular case, the critical 
set Y>a is thus the collection of all line segments joining pairs of eigenvalues of A, and the density /a 
is polynomial of degree n — 3 in each connected component of C \ T,a- In the generic situation where 
no straight line contains more than two eigenvalues of A, one can show that /a is continuous together 
with its derivatives up to order n — 4 (and is discontinuous on dW(A) if n = 3.) 

4. (The nonnormal case) Finally, we consider the most interesting situation where the matrix A G 
M„(C) is not normal. In that case, there is no explicit formula for the numerical density, but the 
problem can be reduced in some sense to the Hermitian case by the following simple observation. For 
any 9 G S 1 = R/(27rZ), let H(9) be the Hermitian matrix defined by 

H(G) = \(e~ i9 A + e i0 A*) = A x cos(fl) + A 2 sin(0) , (2) 

where A x = (A + A*)/2 and A 2 = (A-A*)/(2i). Then Re {e~ ie (Ax, x)) = (H(9)x,x) for all x G dM n . 
Now, if the random variable X is uniformly distributed on <9B n , the distribution of (H(9)X,X) is by 
definition the numerical measure Hh{0) , whereas the distribution of Re (e~ ld (AX, X) ) is easily identified 
as the two-dimensional Radon transform of the numerical measure [ia, evaluated at 9 € [0, 2ir]. We 
thus have 

Kn A {0) = n H{e) , ees 1 , (3) 

where 1Z denotes the two-dimensional Radon transformation. Since the numerical density of H(9) 
is known to be the £?-spline based on the eigenvalues X±(9) , . . . , X n (9) of H(9), we can reconstruct 
the numerical measure ha by inverting the Radon transformation in ([3]), using the the well-known 
backprojection method which plays an important role in tomography [12] . This provides a useful repre- 
sentation formula for the numerical density, as well as an efficient algorithm for numerical calculations, 
see Section H] for more details. 

It is worth mentioning here that the critical set can be conveniently characterized using the 
family of Hermitian matrices H(9) associated with A. Indeed, if we define the eigenvalues Xj(9) in 
such a way that they depend analytically on 9, one can shown that the algebraic curve Ca which 
generates XU is given by 

C A = {e ie (X j (9) + iX' j (e))\je{l,...,n}, 9g[0,tt]} , 

see [261 HU [T7] an d Section [5] below. In the generic case where X\(9) < X2(9) < ■ ■ ■ < X n (9) for all 
9 G [0, 7r], the bitangent set C' A is empty and T,a = Ca- 

As was already mentioned, the numerical density /a is smooth (in fact, real-analytic) on each 
connected component of C\S y 4. The regularity across XU is more difficult to study, but we shall show 
in Section 16.21 that Ja is everywhere of class C n_3 if n > 3 and A G M n (C) satisfies some generic 
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Figure 1: The numerical density Ja is represented for a typical matrix A € M^R), given by ([55]) below. In 
the contour plot (left), the exterior ovate curve is the boundary of the numerical range W(A), and the other 
component of the of the critical set T,a is the interior cuspidal triangle. The three-dimensional plot (right) 
confirms that the numerical density is continuous, positive inside the numerical range, and constant over the 
cuspidal triangle, in agreement with the results of Section [6] 

hypotheses, which exclude in particular the case of normal matrices. In addition, for an arbitrary 
matrix of size n, we shall prove that all derivatives of /a of order n — 2 vanish identically in some 
distinguished regions, which have the following geometric characterization. For any z G C \ T*a, let 
N(z) be the number of straight lines containing z which are tangent to the curve Ca, see ()39p below 
for a precise definition where possible multiplicities are taken into account. It is easy to verify that 
N(z) is constant in each connected component of C \ Ha,, and that N(z) < n [26]. The distinguished 
regions where Ja is polynomial of degree n — 3 (if n > 3) or /a = (if n = 2) are exactly those 
connected components of C \ on which N(z) takes its maximal value n. This remarkable property 
of the numerical density, which is one of our main results, will be established in Section 16.11 The 
geometric condition N(z) = n is always satisfied in the complement of the numerical range, where /a 
vanishes identically, but for many matrices of size n > 3 is it also met in some regions inside W(-A). 
For instance, in the three-dimensional example represented in Fig.[TJ it is easy to verify that N(z) = 3 
if z is outside W(.A) or inside the cuspidal triangle, and N(z) = 1 in the intermediate region where 
the numerical density is not constant. 

At this point, it is necessary to make a connection with the theory of lacunas of symmetric 
hyperbolic systems of partial differential equations in 2 + 1 variables. Given A £ M n (C), we consider 
the following system of linear PDE's in Kt x 1^: 



where A\, A% are as in ([2]) and u = {u\, . . . , u n ) T : 1 x R 2 -> W 1 . The fundamental solution of @ is 
the unique (matrix-valued) distribution E supported in the half-space IR+ x R 2 which satisfies 



One can show that E(t, x) is homogeneous of degree —2 in t and x, it is thus sufficient to consider the 
time-one trace E* = E(l, •), which is a distribution on R 2 . Due to the finite speed of propagation, it 
is well-known that E* is zero outside a compact set of R 2 , but it may also happen that E* vanishes 
identically in some regions inside the domain of influence of the origin. Such regions are called lacunas 
of the hyperbolic system @. 



d t u + Aid xi u + A 2 d X2 u = , 



(4) 



d t E + A^E + A 2 d X2 E = I n 5 t=0 ® ^ =0 . 



(5) 
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The properties of the fundamental solution of symmetric hyperbolic systems have been studied by 
many authors, see e.g. [22J US SI E] . In the particular case of system (j3J), J. Bazer and D. Yen 
have shown that, if one identifies C with M 2 , the singular support of the distribution is contained 
in the critical set T,a, and the (stable) lacunas of system are exactly the regions described above 
where the numerical density Ja is polynomial of degree n — 3. This remarkable coincidence is of 
course not fortuitous. In Section [6.31 we explain it by showing that the fundamental solution E can 
be expressed as a linear combination of derivatives of order n — 1 of a homogeneous extension of 
the numerical density /a- This connection allows us to recover some of the main results of [4], and 
therefore confirms that the numerical measure is a natural quantity attached to the matrix A. One 
might even argue that fiA contains more information than E*, since for instance supp(/iA) = W(A) 
while supp(£ l #) is in general strictly smaller and not necessarily convex, see Sectional Similarly, we 
believe that singsupp(^) always coincide with S^, while singsupp^*) is usually smaller. 

A final question that is worth investigating is the behavior of the numerical measure \ia when the 
size of the matrix A goes to infnity. Here of course, specific assumptions have to be made in order to 
obtain convergence results. Suppose for instance that {A n } n >i is a sequence of complex matrices with 
A n £ M n (C), Tr(^4 n ) = 0, and ||A n || < C for all n > 1. If, for each n > 1, X n is a random variable 
that is uniformly distributed on <9B™, we show in Section [8] that the complex variable (A n X n , X n ) 
converges almost surely to zero as n — > oo. This is reminiscent of the strong law of large numbers in 
probability theory. Under slightly stronger assumptions, we also establish the analog of the central 
limit theorem in this context. Our convergence results mean that [ia is very close to 5 Z when dim(A) 
is large, where z is the barycenter of cr(A). This explains why plotting (Ax,x) for randomly chosen 
points x £ dM n is a very unefficient algorithm for determining the numerical range TU(A) if A is a 
large matrix! 

The rest of the paper is organized as follows. In Section [2j we establish some general properties 
of the numerical measure. Section [3] is devoted to the particular situations where the matrix A is 
Hermitian or normal. The nonnormal case is treated in Sections SHU which constitute the core of the 
paper. In Section 01 we derive a representation formula for the numerical density using the inversion 
of the Radon transformation. Section [5] collects a few results on the geometry of the critical set Y*a, 
which are mainly borrowed from |19[ [26l [T7] . These informations are used in Section [6] to derive an 
explicit formula for the derivatives of order n — 2 of the numerical density, which allows us to obtain 
generic regularity results and to express the fundamental solution of the hyperbolic system (JH) in 
terms of derivatives of the numerical density. To illustrate our results, a few explicit examples are 
treated in Section [71 Finally, we investigate in Section [8] the concentration properties of the numerical 
density for large matrices, and we discuss in Section [9] a possible extension of our results to hyperbolic 
polynomials with an arbitrary number of variables. 

Acknowledgements. This work has benefited of stimulating discussions with several of our col- 
leagues, including Y. Colin de Verdiere, F. Faure, and A. Joye. 

2 General properties of the numerical measure 

In this section, we establish a few general properties of the numerical measure of a complex matrix. In 
particular, we show that [ia is absolutely continuous with respect to the Lebesgue measure on VF(^4), 
and we prove a direct sum formula which will be useful later. 
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2.1 Support and regularity properties 

We first show that the support of the numerical measure always coincides with the numerical range 
of the matrix. 

Lemma 2.1 For any A G M n (C), one has supp(//A) = W(A). 

Proof. If V = C \ W(A), then <S> A \V) = 0, hence n A (V) = a(<Z> A \V)) = 0. This shows that 
supp(/iA) C W(A). Conversely, if V C C is any open set such that V D W(A) / 0, then & A l (V) is a 
nonempty open subset of dM n , hence ha(V) = a^^iV)) > 0. Thus W(A) C suppi^). ■ 

Our next goal is to locate the singular support of [ia- We recall that x G dM n is a regular point 
of &a if the differential map d x &A '■ T*dM n — > C is onto. Otherwise, we say that x is a critical point. 
The following characterization will be useful: 

Lemma 2.2 pEH [[7] Lei A G M n (C). 

i j j4 point x G <9B n is a critical point of the numerical map &a if o-nd only if x is an eigenvector of 
the Hermitian matrix H(6) defined in ([2]) for some 8 G [0, 7r]. 

The differential of &a vanishes at x G <9B n i/ and only if x is an eigenvector of both A and A* . 

In other words, the range of the differential d x &A has (real) dimension 1 if and only if x is an 
eigenvector of H{0) for a unique 9 G [0,7r), and is reduced to {0} if and only if x is an eigenvector of 
H(9) for all 9 G [0,7r]. The proof is neither new nor difficult, but we shall repeat it here in order to 
introduce some notation that will be needed later on. 

Proof. Since ^^(e^x) = &a(%) for all 9 G [0,2-7r], we can consider the numerical map as acting on 
the quotient space dW l / S l ~ CP n_1 [16]. Thus, to detect the critical points of we study the 
reduced map $a ■ dM n /S 1 -)• C defined by 

$a(W) = $a{x) = (Axx, x) + i(A 2 x, x) , x G <9B n , 

where [x] = {e ie x \ 9 G S 1 } and A\,A 2 are the Hermitian matrices introduced in ([2]). 

If x G dM n , the tangent space to dM n /S 1 at [x] is just the (2n — 2)-dimensional afhne subspace 
{[x + y] | y G C n , (x,y) = 0}. Thus, using the definition above of $a, it is straightforward to verify 
that, for all y G C n with (x,y) = 0, one has 

^d [x] $ A (y) = (A lX \y) + i(A 2 x\y) = {v x {x)\y) + i{v 2 {x)\y) , (6) 

where (x\y) = Ke(x,y) = (Rex)*(Rey) + (Imx)*(Imy) denotes the real scalar product in C n ~ IR 2 ™, 
and 

v i(x) = A\x — {A\x, x)x , v 2 (x) = A 2 x — (A 2 x, x)x . 

Of course, replacing A\x, A 2 x with t>i(x), v 2 {x) has no effect in © since (x,y) = 0, but after this 
substitution we can let y run over the whole of C n without increasing the range. So our task is reduced 
to computing the rank of the IR-linear map y i— > (vi(x)\y) + i(v 2 (x)\y), which is just the rank of the 
2x2 matrix 

r>/Vi - f (Mx)\vi{x)) {vi{x)\v 2 (x))\ 
{ ' \{v 2 {x)\ Vl {x)) (v 2 (x)\v 2 (x))J ■ [t) 
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By the Cauchy-Schwarz inequality, the positive matrix D(x) is singular if and only if there exists 
6 G [0, 7r] such that v\{x) cos 9 + V2{x) sin 6 = 0, which exactly means that x is an eigenvector of H(9). 
Moreover, D(x) = if and only if v\{x) = V2(x) = 0, which is equivalent to saying that x is an 
eigenvector of both A\ and A2 , hence of both A and A* . m 

Let T*a C C denote the set of all critical values of &a, namely = <I>^4(r(^4)) where T(A) C dM n 
is the set of all critical points of <&a- Our next result is: 

Lemma 2.3 If A £ M„(C) ; then sing supp(/u y i) C S^. 

Proof. If the numerical range ^(^4) is reduced to a line segment or to a single point, then Y>a = W(A), 
hence singsupp(^ y i) C supp(A) = T,a by Lemma l2.1i Thus, we assume from now on that VF(>1) has 
nonempty interior. By Sard's lemma, the critical set T,a is then a compact subset of with zero 

Lebesgue measure. We have to show that there exists a smooth density function Ja > such that 
d/iA(z) = Ja(z) dz on C \ XU- Clearly, we must have Ja = on C \ W(A). 

If z e W{A) \ S A , then M z := ^{z) is a compact submanifold of dM n of codimension 2, which 
depends smoothly on z. Using classical arguments, involving a partition of unity and the Implicit 
Function Theorem, it is not difficult to verify that, for any continuous function <j) with supp((/>) C 
W{A) \ T, A , one has 

/ 0((Ax,x))da(x) = — I <t>(z){! dz , 

where uj n = 2ir n /((n — 1)!) is the total measure of dM n , v is the {2n — 3)-dimensional Euclidean measure 
on the submanifold M z , and A(x) = det-D(x) where D(x) is the 2x2 matrix defined in (|7|). Remark 
that 2A(x) 1 / 2 = Xi(x)\2(x), where Ai(x), A2(x) are the singular values of the differential map d x &A- 
In view of ([1]), we conclude that dfiA(z) = /a(z) dz on C \ S^, where 

fA(z) = -i- / z € C \ Tia ■ (8) 

It is easily verified that the density /a is smooth and strictly positive on W(A) \ Y*a- ■ 

The results obtained so far are summarized in the following proposition, which also asserts that 
the numerical measure is absolutely continuous with respect to Lebesgue's measure on iy(j4). 

Proposition 2.4 Let A £ M„(C). 

1) If the numerical range W{A) has nonempty interior, the numerical measure ha is absolutely con- 
tinuous with respect to the (two-dimensional) Lebesgue measure on C ~ M. 2 . The numerical density 
fA = d^A/dz is smooth outside the critical set T, A - 

2) If A is a nonscalar Hermitian matrix, then W(A) C M and the numerical measure is absolutely 
continuous with respect to the (one-dimensional) Lebesgue measure on M. The numerical density 
fA = dfiA/dx is smooth outside the spectrum a (A). 

Remark 2.5 As is explained in the introduction, Proposition \2.J\ covers all interesting cases. Indeed, 
if the numerical range has empty interior, then either W (A) is reduced to a single point, in which case 
A is a scalar matrix and [ia is just a Dirac mass, or W (A) is a line segment of nonzero length, in 
which case A can be reduced to a nonscalar Hermitian matrix by a simple affine transformation. 
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Proof. Using the same notations as in Lemmas 1 2 . 2 1 and [2T3| we observe that F(A) = {x G dM n \ A(x) = 
0}, where A(x) = det-D(x) is a polynomial in the 2n variables Rexj,Imxj (i = 1, . . . , n). Thus one of 
the following two situations must occur: 

1) T(A) is an algebraic submanifold of dM n of codimension at least 1. By Sard's lemma, this is the 
case if and only if VF(j4) has nonempty interior. In that situation, since we already know that fiA has 
a smooth density outside the critical set XU, we only need to show that /m(Xa) = 0. Given e > 0, 
let T e (A) = {x G dM n | dist(x, r(^4)) < e}, where "dist" denotes here the geodesic distance on the unit 
sphere. We decompose 

$7 X (E A ) = f*7 1 (SA)nr e (A) > )u('*7 1 (E il )nr e (A) c> ) = E 1 (e)UE 2 (e) , 



where T e {A) c = dM n \ T e {A). Since a(T(A)) = 0, we have a(Ei(e)) < a(T e (A)) -> as e -> 0. 
Moreover, the proof of Lemma 12.31 shows £>2( e ) is a codimension two submanifold of <9B n , so that 
a(E2(e)) = for any e > 0. Using the definition of the numerical measure, we conclude that /i^(S^) = 
a(^ 1 (SA)) = 0. 

2) T(A) = dM n . This is the case if and only if W(>1) has empty interior, and without loss of generality 
we can then assume that the matrix A is Hermitian. Since W(.A) C R, it is more natural here to 
consider $^ as a map from <9B n into R. If we do that, then repeating the proofs of Lemmas 12.21 and 12. 31 
we easily find that the critical points of &a are exactly the eigenvectors of A. Moreover, the numerical 
measure has a smooth density on R\ cr(A), and is absolutely continuous with respect to the Lebesgue 
measure dx if A is not a scalar matrix. We skip the details here, because the Hermitian case will be 
treated in full details in Section [3] below. ■ 

To conclude this section, we show that the numerical density is strictly positive in the interior of 
iy(j4). A little care is needed in the formulation of that result, because as we shall see in Section [4.21 
the numerical density need not be a continuous function. 

Proposition 2.6 If zq G C is an interior point ofW(A), then 

liminf \ ^ A ({z G CI \z - z \ < e}) > . (9) 

Proof. If zq is an interior point of iy(^4), it is shown in |1T|, Proposition 2.11] that the preimage 
$^ 1 (zo) contains at least one regular point xq. Let V be an open geodesic ball centerd at xq G dW n 
whose closure does not intersect T(A). Proceeding as in the proof of Lemma 12.31 we find 

H A {{zeC\\z-zo\ <e}) > — / (/ - y 1 dz . 



Wn J\z-z \<e UM z nv 2A(x)V2 

If e > is sufficiently small, the integral inside the curly brackets is a smooth and positive function of 
z, and ([9]) follows. ■ 

2.2 The direct sum formula 

Let p, q G N* and n = p + q. Given A G M p (C) and B G M g (C), the direct orthogonal sum of A and 
B is the matrix A B G M n (C) defined by 

a n /AO 

A(BB ={o B 

In this situation, we have a formula for the numerical measure ha®b in terms of ha and [is- 
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Proposition 2.7 For any <fi £ C°(C), we have 

>(z)dfx AeB (z) = f f f 1 ^tz' + (l-t)z")tP- 1 (l-ty- 1 dtdfi A (z')dv B (z") , (10) 

B{p,Q) Jc Jc Jo 

where B(p,q) is Euler's beta function 

Proof. Any unit vector x £ dM n can be written as 

x=( ^ U ) 

where u £ dW, v G SB 9 , and t £ [0,1]. Up to negligible sets, the map (u,v,t) t-t x defines a 
diffeomorphism from dM p x dM q x [0, 1] onto dM n . With this parametrization it is not difficult to verify 
that the Euclidean measure on dM n has the following expression 

da n (x) = ^tP-^l-ty^dtdapi^daqiv) . 

Equivalently, since o~ n (dM n ) = uj n = 2-7r n /T(n), the normalized Euclidean measure satisfies 

d ^ x) = -^J ) tP ' 1 ( 1 - t ) q ' ldtd ^p( u ) d ^( v ) • 
Thus, using definition (pQ) and the fact that {(A © B)x, x) = t(Au, u) + (1 — t)(Bv, v), we easily obtain 



(z)dHA(BB(z) 



f I I 4>(t(Au,u) + (l-t)(Bv,v))t p - 1 (l-t) q - 1 dtda q (v)da p (i 

^ // I' 4>(tz' + (l-t)z")t p - 1 (l-ty- 1 dtd f x A (z')d f i B (z") , 
3, q) Jc Jc Jo 



B(p 

which is the desired result. ■ 

As an application, if we choose <j>(z) = e~^' z in Proposition 12. 7| we obtain the following relation 
between the Fourier transforms of the measures /j,a, fJ-B and ^a®b '■ 

A*bb(0 = -^7^ / MtO Ab((i - 1)0 t p -\i - tf- 1 dt, e e M 2 . (li) 

o(Pi Q) Jo 

The formula given in Proposition 12.71 can be generalized in a straightforward way to a direct sum 
with an arbitrary number of terms. Assume that A = A\ © • • • © Ak where Aj £ M p . (C), so that 
A £ M n (C) with n = pi + ■ ■ ■ + pt- Setting p = {pi,... ,Pk), we denote 



B(p) 



/ tr 1 ---tr i dt 1 ...dt^ 1 = ^ 

Jd^ l (Pi 



r(p i) • • • r( Pk ) 
+ -"+Pk) 
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where t k = 1 — (t\ + ■ ■ ■ + t k -i) and D k _\ denotes the (k— l)-dimensional simplex 

D k -x = {(*!,..., i fc _i) e M^T 1 1 H h < l} • (12) 

Using (jlOp and proceeding by induction over k, we easily obtain the general formula 

/ <j>{z) dflAiz) = I [ cfrihZ! + ■■■+ t k z k ) if' 1 

dti . . . dt fc -i d/UAi^i) • • • dfi Ak (z k ) , (13) 
where it is understood again that t k = 1 — (ti H + ifc-i). 

3 The numerical density of a normal matrix 

If A € M„(C) is a normal matrix, the numerical measure [ia is entirely determined by the spec- 
trum cr(A) = {Ai,...,A n }. Indeed, we know that A is unitarily equivalent to the diagonal matrix 
diag(Ai, . . . , A n ), and that a unitary conjugation does not affect the numerical measure. Using this 
observation and the direct sum formula of Section 12.21 we shall prove that the numerical density of A 
is a piecewise polynomial function, which can be characterized as a multivariate S-spline whose knots 
are the eigenvalues of A. We begin with the important particular case where all eigenvalues of A are 
colinear. 

3.1 The Hermitian case 

If A G M n (C) is a Hermitian matrix, then IU(yl) C M and the numerical measure [la is therefore 
supported on the real axis. Assuming that A is not a multiple of the identity matrix, we show in 
this section that fiA is absolutely continuous with respect to Lebesgue's measure on R, and we give a 
simple characterization of the numerical density f A = d/x^/dx. The result is: 

Proposition 3.1 If A £ M n (C) is a nonscalar Hermitian matrix, the numerical density fA '■ K — > M. + 
is the normalized B-spline of degree n — 2 whose knots are the eigenvalues of A. 

To make the statement clear, we briefly recall the definition and some elementary properties of 
the classical l?-splines [7]. If X%, . . . , X n 6 R are pairwise distinct, the (n— l) th divided difference of a 
continuous function g at the points Ai, . . . , A n is the quantity 

^-^--ird^y (14) 

This is the leading coefficient of the unique polynomial of degree at most n— 1 which agrees with g at 
the points Ai, . . . , A n . It is easy to verify that the right-hand side of (|14p is a completely symmetric 
function of the variables Aj. If g £ C n-1 (1R), the divided difference can be extended by continuity to 
arbitrary (not necessarily distinct) values of Ai, . . . , A n , and we have the integral formula: 

<J n - 1 [Ai,...,A n ] 5 = f g( n - 1 Hh\i + --- + t n \ n )dt 1 ... dt n _! , (15) 
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where D n -\ is the (n— l)-dimensional simplex defined in (|12|) and t n = 1 — (t\ + • • • + in-i)- In what 
follows, we shall always assume that the set S = {Ai, . . . , A n } is not reduced to a single point, so that 
the (n— l) th divided difference is well-defined as soon as g is of class C n ~ 2 in a neighborhood of S. 

With these notations, the normalized B -spline of degree n — 2 with knots Ai, . . . , A n is the function 
B : R -»• R defined by 

= i?[Ai, . . . , A n ](x) = (n — 1) <5 n-1 [Ai, . . . , A n ](- - x)™~ 2 , ieK, (16) 

where (• — x)+ -2 denotes the map y h-> max(0,y — x) n_2 . If Ai < • • • < A n , it is not difficult to show 
that B(x) vanishes identically outside [Ai, A n ], and coincides with a polynomial of degree at most n — 2 
on each nonempty interval (Xj, Aj+i). Moreover, if rrij denotes the multiplicity of Xj in S, one can 
verify that B(x) is continuous at x = Xj together with its derivatives up to order dj = n — 2 — rrij, 
provided dj > 0. If dj = —1, then B{x) is discontinuous at Aj. Finally, we shall see below that B(x) 
is positive on (Ai, A n ) and that L B(x) dx = 1. 

Proof of Proposition 13.11 Let A G M n (C) be a normal matrix with eigenvalues X\, . . . ,X n . To 
compute the numerical density, we can assume without loss of generality that A is diagonal, namely 
A = Ax © • • • © A n with Aj = Xj G Mi (C) . Thus we can use the direct sum formula (I13D with k = n 
and p = (1, . . . , 1). Since B{p)~ l = (n — 1)! and ^a- = 8\. for j = 1, . . . , n, we obtain the relation 

I <t>(z) dfi A (z) = (n - 1)! / 0(tiAi + • • • + t„A n ) dti . . . dt n _! , (17) 

for any continuous function <p : C — > C. 

Assume now that ^4 is Hermitian, so that Ai, . . . , A n 6K, and that o~(A) is not reduced to a single 
point. If B = JB[Ai, . . . A n ] is the normalized i?-spline defined by (|16H . we claim that 

/ <^(x)5[Ai, . . . , A n ](x) dx = (n - 1)! / <^(tiAi + • • • + t n A n ) dti . . . dt n _i . (18) 

Indeed, it is clearly sufficient to prove (fT8j) for compactly supported functions 4> G C°(1R). Moreover, 
since both members of (fT8|) depend continuously on Ai, . . . , A n , we can also assume that the eigenvalues 
of A are all distinct. In that case, it follows immediately from (I14D . (I16p that 

/ 0(a:)B[A ll ... J A n ](x)da: = (n - 1) d^^Ai, . . . , A n ]$ , 

where $(y) = L <j)(x)(y — x)^T 2 dx. Since $( n_1 ) = (n — 2)! cf>, we can use the integral formula (fT5|) to 
evaluate the divided difference in the right-hand side, and we obtain (I18p . 

Now, comparing (fT7|) and ([15]) , we conclude that is absolutely continuous with respect to 
Lebesgue's measure on R, and that the numerical density /a = dfiA/dx is precisely the normalized 
B-spline B[Xi, . . . , A ra ]. Incidentally, the argument above shows that B is positive on its support (see 
Proposition 12. 6p and that f R B(x)dx =1. ■ 

3.2 The quasi-Hermitian case 

We say that a matrix A G M n (C) is quasi-Hermitian if the numerical range VK(^4) C C has empty 
interior, i.e. W(.A) is a single point or a line segment. In such a case, there exist z G C, 9 G S , 
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and a Hermitian matrix H such that A = zl n + e lS H (in particular, A is normal). Indeed, if we 
choose z,6 such that W(A) C z + e* 9 R, the matrix H = e~ i6 (A - zl n ) satisfies W(H) C R and is 
therefore Hermitian. Since (Ax,x) = z + e t9 {Hx,x) for all x E <9B n , it is clear from ([1]) that the 
numerical measure \ia is just the image of [in under the affine isometry w i— )• z + e l6l u7. Combining 
this remark with Proposition 13.11 we thus obtain a precise characterization of the numerical measure 
of any quasi-Hermitian matrix. 

3.3 The normal case 

Finally, we consider the case of a normal matrix A £ M n (C) whose numerical range W(-A) has 
nonempty interior. This, of course, is possible only if n > 3. By Proposition 12.41 the numerical 
measure \xa is absolutely continuous with respect to Lebesgue's measure, and we have the following 
characterization of the numerical density /a = d/^/dz: 

Proposition 3.2 If A € M n (C) is a normal matrix whose numerical range W(A) C C has nonempty 
interior, the numerical density /a ■ C — > R+ is the bivariate B-spline of degree n — 3 whose knots are 
the eigenvalues of A. 



The reader is referred here to the work of W. Dahmen [6] , where multivariate i?-splines are defined 
and studied in detail. To make the connection with the numerical density of a normal matrix, we use 
the relation (|17p . which corresponds to formula (2.2) in [6]. In the rest of this section, we assume that 
A £ M„(C) is a normal matrix whose eigenvalues Ai, . . . ,X n are not colinear, and we often identify 



Proof of Proposition 13.21 Assume first that n = 3. Then W(A) C C ~ R 2 is the 2-simplex with 
vertices Ai, A2, A3, and using the change of variables w = tiAi + ^A2 + (1 — t\ — £2^3 in (fT7|) we easily 
obtain 



This shows that the numerical measure \ia is uniformly distributed on W(A). The numerical density 
is thus a multiple of the characteristic function of W(.A), which (by definition) is the bivariate -B-spline 
of degree zero with knots Ai, A2, A3. 

We now assume that n > 4. Then we can choose n vectors v\,...,v n G M n_1 such that 

1) The (n— l)-simplex S C R"" 1 with vertices v\, . . . , v n has unit volume; 

2) Pvi = Aj for i = 1, . . . , n, where P : R n_1 — >• R 2 is defined by P{x\, . . . , x n -\) = (x\,X2). 

This (elementary) claim is proved in [6l Section 2]. The simplex S C R n_1 is not uniquely defined, 
but any choice satisfies PS = W{A) = conv(<r(^4)). Returning to (fT71) . we have 



where the second equality is obtained by applying the change of variables (t%, . . . , i n -i) w = 
t\V\ + • • • + t n v n , with t n = 1 — (ti + • • • + t n -i). This shows that the numerical measure is the image 
under the projection P of the Lebesgue measure on the simplex S C R n_1 . Given z S C ~ R 2 , the 
numerical density /a(-z) is thus the (n— 3)-dimensional measure of the simplex P~~ 1 z n S. This is 
precisely the definition given in [6] of a bivariate i?-spline of deg ree n — 3 with knots Ai , . . . , A n . ■ 



the complex plane C with R 2 . 





' 4>{Pw) dw 
s 
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We conclude by listing a few properties of the numerical density Ja, which follow from [6 3 The- 
orem 4.1]. Let Ha C C be the union of all line segments joining pairs of eigenvalues of A. Using 
Lemma 12.21 it is straightforward to verify that T*a is exactly the set of critical values of the numerical 
map &a- Then /a is a polynomial of total degree n — 3 in each connected component of C \ S^. 
Moreover, if £ C C is a straight line passing through a pair of eigenvalues of A, the numerical density 
is continuous across £ n W(A) together with its derivatives up to order n — 2 — m > 0, where m > 2 is 
the number of eigenvalues of A (counted with multiplicities) which belong to £. If m = n — 1, then /a 
is discontinuous on t n W(-A). In particular, in the generic case where no straight line contains more 
than two eigenvalues of A, the numerical density Ja is of class C" 1-4 if n > 4. 

Remark. In the Hermitian and normal cases, the fact that the numerical density /a is a projection 
of the characteristic function of a convex set, which is a log-concave function, implies that the density 
is itself a log-concave function. We thus have the following inequality 

f A (Xz + (1 - X)z') > f A (z) x f A {z') l - X , (19) 

for all z,z' G C and all A G (0,1). This follows from the Prekopa-Leindler inequality, see e.g. [TU1 
Section 9]. We warn the reader that this property does not extend to nonnormal matrices, as we can 
see already from the two-dimensional case considered in Section [4.1i 

4 The Radon transform of the numerical measure 

Our purpose in this section is to derive a representation formula for the numerical density of a non- 
normal matrix A G M n (C). Our approach is based on a natural expression of the Radon transform 
of the numerical measure ha in terms of the Hermitian matrices H(9) defined in §2$. By definition, 
the Radon transform of fiA is the family IZfiA = {7^a(#) I £ 5" 1 }, where 1Zha{6) denotes the Borel 
measure on R defined by 

(Kha(0))(I) = li A {{z e C| Re(e- ie z) G /}) , 

for any open set I C R. In other words, 1Z[ia{Q) is the image of the measure [la under the orthogonal 
projection in C ~ M. 2 onto the line e i0 M. The fundamental observation is: 

Proposition 4.1 For any 9 G S 1 = M/(27rZ) ; one has TZ^a{Q) = ^H(6)- 

Proof. The definition ([2]) implies that Re (e~ l9 (Ax, x)) = (H(0)x,x) for any x G dM n . Thus, for any 
open set I C R, we have 

(Kha(0))(I) = a({x G dM n \ Re (e~ ie (Ax , x)) G J}) 

= a({x G 8M n \ (H(0)x,x) G I}) = ^ H (e)(I) , 

which proves the claim. ■ 

Proposition 14.11 shows that the numerical measure of an arbitrary matrix A G M„,(C) is entirely 
determined by the one-dimensional measures associated with the Hermitian matrices {H(0) \9 G S 1 }. 
If we assume that W(A) has nonempty interior, which is always the case if A is nonnormal, the matrix 
H(9) is nonscalar for every 9 G 5 1 and Proposition 13. II show that its numerical density is the S-spline 
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B[Xi(9), . . . , X n (9)], where Xi(9), . ■ . , X n (9) are the eigenvalues of H(0). In that case, the result of 
Proposition 14. II can be stated in the following equivalent form 

/ f A (e i6 (x + iy))dy = B[Xx(9), . . . , X n (9)](x) , (20) 

JR 

where equality holds for all 9 £ S 1 and almost all x G M since the numerical density f a belongs to 

lHc). " 

Our goal is to invert the Radon transform (|20p to obtain a representation formula for the numerical 
density Ja- The general results established in |12j show that 

(Jf A )(x + iy) = -jj- [ B[X 1 (9),...,X n (9)](xcos9 + ysm9)d9 , (21) 
4vr J s i 

where J = (— A) -1 / 2 is the Riesz potential defined by (Jf)(z) = ^ J c \z — z'\~ 1 f(z') dz' . The idea 
is thus to apply the nonlocal operator (— A) 1 / 2 to both sides of (|21|) . but since we are not dealing 
with smooth functions we have to differentiate in the sense of distributions. As a preliminary remark, 
if f(x,y) = g(xcos9 + ysin9) for some test function g : R — > R and some fixed 9 E S , a direct 
calculation shows that (Af)(x,y) = g"(xcos9 + ysin9), and a standard interpolation argument allows 
us to conclude that (— A) 1 / 2 /(x, y) = Tig' (x cos9 + ysin9), where Tig' denotes the Hilbert transform 
of the derivative g' . Using this observation, we easily obtain the representation formula 

f A (x + iy) = ^f HB'[Xi(9),...,X n (9)}(xcos9 + y S in9)d9 , (22) 
4vr J s i 

where both sides define integrable functions of z = x + iy E C, and equality holds almost everywhere. 
So we have shown: 

Proposition 4.2 If the numerical range of a matrix A G M n (C) has nonempty interior, the numerical 
density of A can be represented as in (|22p for almost all (x,y) £ M 2 . 

In the rest of this section, we shall apply Proposition 14.21 to compute the numerical density of a 
two-dimensional nonnormal matrix. We shall also consider the interesting particular situation where 
the numerical density is radially symmetric, in which case the representation formula takes a simpler 
form. Proposition 14.21 will be used again in Section [6] to derive some important properties of the 
numerical measure in the general case. 



4.1 The two-dimensional case 

Let A G M2(C), and assume that the numerical range VK(^4) has nonempty interior. As is well-known 
|15| . VF(^4) is then a filled ellipse, and without loss of generality we can assume that this ellipse is 
centered at the origin and that its major axis is aligned with the real axis of the complex plane. In 
that case, up to a unitary conjugation, the matrix A has the following form 



.4 



-c 26 N 
c 



? ) , W(A) = {x + iy£c\^ + ^<l} , (23) 

where b > 0, c > 0, and a = (b 2 + c 2 ) 1 / 2 . For 9 G S 1 = R/(2ttZ), let 

H{9) = he-»A + e»A*)= H™ 9 ^ . 
w 2 \ be c cos 9 ) 
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The eigenvalues of H(9) are ±A(0), where X(6) = (b 2 + c 2 cos 2 6*) 1 / 2 . It follows that 

B[-\(e),x(e)](s) = ^i hmxe)] ( s ) . 

According to (|22p . we have to differentiate this expression with respect to s, which yields a linear 
combination of Dirac masses, and to apply the Hilbert transformation, which is the convolution with 
the distribution p. v. — . We obtain 

^ TVS 

HB'[-\(0),\(0)](s) = -J—fp.v.^— --p.v. ^r) = -p.v. ' 



2ir\(6) V*""s + A(#) t ""s-X(9)J vr ' X(9) 2 - s 2 ' 
hence 

If 1 

f A (x + iy) = -j-o-p.v. / 7 2 ■ 2 2fl — 7 a - ! ^~7JW d61 - ( 24 ) 

47T 7s 1 " + c cos 9 — [x cos p + y sin vy 

It remains to compute the right-hand side of (|24|) . which is a simple exercise in complex analysis. 
Setting z = x + iy and w = e 2ld , we first observe that 

6 2 + c 2 cos 2 #- (xcos# + ysin#) 2 = — ((c 2 - z 2 ) + 2w{a 2 + b 2 - \z\ 2 ) + w 2 (c 2 - z 2 )) , 

hence 

fA(z) = — P-V. * 71 2 , - p 2 3 T-T2V- 27-3 =2T dw ■ ( 25 ) 

Z7r z yitu|=i ( c ~~ z J + 2w(a z + b z — \z\ z ) + ur(cr — z z ) 



In (I25p . the roots of the denominator are 

w± = -^-^(-a 2 -b 2 + \z\ 2 T vV + & 2 " kl 2 ) 2 " (c 2 " ^ 2 )(c 2 - ^ 



c 2 



1 



c 2 — z 2 



,2 



-a 



6 2 + |z| 2 T V4(a 2 ^ 2 - & 2 ^ 2 - a 2 y 2 ) . (26) 



We can therefore distinguish between two cases: 

1. The point z = x + iy belongs to the interior of W(-A). Then the expression under the square root 
is positive, and it is easy to verify that \w+\ > 1, \w—\ < 1 (in the limiting case where z = ±c is a 
focus of the ellipse, one can set w- = and w+ = oo.) Thus the principal value in (|25p is not needed, 
and the residue theorem shows that 

11 11 
fA{z) = n(a 2 + b^-\z\ 2 ) + (c 2 -z z )w^ = ^vW-& 2 x 2= P? ' (27) 

2. The point z = x + iy lies outside W(-A). Then the expression under the square root in (|26p is 
negative, and one verifies that \w±\ = 1. Thus the integrand in (I25p is holomorphic outside the unit 
circle {\w\ = 1} and decreases like l/|w| 2 at infinity. It follows that 

h ^ := i^ f H=r (c 2 - z 2 ) + 2w{a 2 +b 2 - |z| 2 ) + w 2 (c 2 - z 2 ) dW = ° ' 

for all r + 1, hence f A (z) = \{I Z {1+) + 4(1-)) = 0. 

Summarizing, we have shown that the numerical density of the matrix (|23p is the function f A G 
L X (C) defined by (|27l) inside W(A), and vanishing identically outside W(A). Note that singsupp(/iA) = 
dW(A) = Y>a, in agreement with Lemma [2731. an d that Ja(z) blows up when z converges to the bound- 
ary of VF(^4) from inside. In particular f A is not log-concave, in contrast to what happens when A is 
normal. 
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4.2 The radially symmetric case 



It sometimes happens that the numerical range of a matrix A G M n (C) is a disk in the complex plane 
and that the numerical density is radially symmetric about the center. In such a case, it is possible 
to obtain a representation formula which is simpler than (|22p . Indeed, assume that the disk W(>1) is 
centered at the origin, and let R > denote the numerical radius of A (or any larger positive number). 
We set 

f A (z) = F A (R 2 - \z\ 2 ) , zeC, \z\<R, (28) 

where Fa ■ [0, R?] — > 1R+ has to be determined. Since the numerical measure is invariant under 
rotations about the origin, it follows from Proposition 14.11 that the projected measure ^H(e) does not 
depend on 6. In analogy with (|28p . if fn denotes the numerical density of H(8) for any 6 £ S 1 , we set 

f H (x) = F H {R 2 - x 2 ) , x G M , \x\<R. (29) 

We then have the following result: 

Proposition 4.3 If the numerical density of a nonscalar matrix A G M n (C) is radially symmetric, 
the functions Fa, Fh defined in (j28j) . ([29]) satisfy the relations 

F H (t) 



[ F A (t-s)-^ds, F A (s) = -^-[ F H (s-t)\dt, (30) 
Jo V s n ds Jo Vt 



fort,s G [0,R 2 ]. 



Proof. By the definition of the Radon transformation, we have fij{x) = L fA(x + iy) dy for all x G 
Using (|2"gj) . ([2"9"]) and the support property, we thus find 

F H (R 2 -x 2 )=j F A {R 2 -x 2 -y 2 )dy = / F A (R 2 - x 2 - s)-j= ds , 



2 <R 2 -:r 2 JO 



for all x G [— R, R]. Setting t = R 2 — x 2 , we obtain the first relation in (|30p. So far, we have shown 
that Fh = ^^IFa, where I is the Riesz potential 



(//)(*) = 4= f f(t-s)^=ds , t>0 
V 71 " Jo Vs 



Now, it is well known that (I 2 f)(t) = Lf(s)ds, see e.g. [121 Chapter V.5], thus the second relation 
in (1301) follows from the first one. ■ 



Examples. Let A = (aij) be a complex matrix, and assume that there exists a nonzero integer k 
such that aij = whenever j — i ^ k. Then iy(j4) is a disk centered at the origin, and the numerical 
density is radially symmetric. Indeed, given any 6 G S 1 , the map 





M 




( e ie xx \ 


X = 




!->• Xg = 










\e in6 x n j 



is a measure-preserving isomorphism of the unit sphere dM n , and our assumption on A implies that 
{Axq,xq) = e lk9 {Ax,x) for all x G dM n . This proves that a rotation of angle k6 about the origin 
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does not affect the numerical measure (jla- Since 6 is arbitrary, the numerical density /a is necessarily 
radially symmetric. 

The simplest example in this category is the 2x2 Jordan block 

m = (3D 



K °, 

Here W(A 2 ) = {z G C | \z\ < 1}, and applying (J27|) with a = b = 1 we find 

Mz) = ^7r=pF 1{W<1> ' (32) 

Alternatively, since the eigenvalues of H 2 = \{A 2 + A£) are ±1, we have fjj 2 = \ l[-i,i] and applying 
Proposition 14.31 we easily obtain 



As a more interesting application, consider the 3x3 matrix 

/0 a 0\ 

A 3 = b , (33) 




where a, b E C and |a| + \b\ > 0. Multiplying A3 with a positive constant, we can assume that 
\a\ 2 + \b\ 2 = 4. Then W{A 3 ) = {z G C | \z\ < 1}, and 



1 1 + y/1- |z| 2 
/a 3 (^) = - log n 1 {0<[«[<l} • ( 34 ) 

7T |Z| 

Indeed, the eigenvalues of ^3 = ^(^3 + A%) are and ±1, hence fH 3 {x) = B[—l, 0, l](x) = (1— 
Applying Proposition 14.31 we obtain (|34p by a straightforward calculation. Alternatively, if we take 
(a, 6) = (0,2), we can use the direct sum formula (fTUj) to deduce (j34l) from (|32j) . Here again, we 
have singsupp(/iA 3 ) = cW(^3) U {0} = Sa 3 , in agreement with Lemma [2731 Remark that fA 3 (z) = 
C((l— |z|) 1//2 ) as \z\ — > 1—, so that the singularity of the numerical measure at the boundary is weaker 
than it was for A 2 . This reflects the fact that fH 3 (%) = (1 — \ x \)+ is Lipschitz continuous, whereas 
fH 2 i x ) = l[-i,+i]( x ) had jump discontinuities. However, we observe that the (logarithmic) singularity 
of /j4 3 at the origin is much stronger than the (square root) singularity at the boundary. As we shall 
see in Section [5j this is a nongeneric concentration phenomenon due to the fact that the component 
of the critical set associated with the eigenvalue of H(8) is reduced to a single point. We have 
here the rare instance of an unbounded numerical density for a matrix of size n > 3. 

Another interesting conclusion that can be drawn from this example is that the numerical density 
HA does not determine the matrix A £ M n (C) up to unitary conjugations if n > 3. Indeed, if we set 
(a, b) = (0,2) and (a, b) = (y/2, \/2) in the definition (|33l) of .A3, the resulting matrices are not even 
similar, yet they have the same numerical density, given by (pE 



5 The geometry of the singular set 

We know from Lemma 12.31 that the numerical density of a matrix A G M n (C) is smooth outside the 
set T,a of all critical values of the numerical map &a- I n this section, we describe a few geometrical 
properties of the singular set 'Ea which will be needed to formulate our main results in Sectional We do 
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not claim much originality here: the material of this section is essentially borrowed from [19 1 HI [16 ] 117] . 
and is reproduced below for the reader's convenience. 

As was shown by Kippenhahn [19], the singular set T,a has a natural description in terms of the 
eigenvalues Xi(6), ■ ■ ■ , X n (9) of the Hermitian matrices H{9) defined in (J2j) . To see that, we first recall 
that these eigenvalues can be numbered in such a way that they are real-analytic functions of 9 G R, 
see [26j. By analyticity, any two eigenvalues either coincide for all 9 G R or cross at most a finite 
number of times on each compact interval. As a consequence, there exists an integer m <n such that 
H(9) has exactly m distinct eigenvalues for all £ [0, ir) \ 0, where C [0, tt) is a finite set. Since 
H{9 + 7r) = —H(9), it follows that we can number the eigenvalues in such a way that 

a(H(9)) = {Ai(0),...,A m (0)} , for all 9 G R , 

where X\(9), . . . ,X m (9) are pairwise distinct and have constant multiplicities outside the crossing set 
+ irZ. Let r : {1, . . . , m} — > {1, . . . , m} be the permutation defined by 

{A T(1) (0),...,A r(m) (0)} = {^(9 + *),... ,-^(9 + , (35) 

for any 9 G [0, 7r) \ (hence for all G R). If we decompose r into disjoint cycles ^i, . . . we 
can associate to each cycle its length £j and its multiplicity mj, the latter being defined as the 
multiplicity of Xj(9) as an eigenvalue of H{9) for any j G ^j. By construction, we have 

£\ + ■ ■ ■ + £k = m > an d 4mi + • • • + ikfnk = ft • 

Moreover, if j G ^j, then Xj(9 + £ jit) = {— Xj(9), hence Xj(9) is periodic with period £jtt if £j is 
even and 2£jir if £j is odd. Note however that these periods are not necessarily minimal. 

Now, we associate to each cycle of the permutation r a closed curve Cj C C defined by 

Cj = {e l0 (Xj(9) + iX'^9)) | G R} , (36) 

where j G {1, . . . , m} is any element of the cycle ^j. Since Xj{9 + 2£jir) = Xj(9), it is clear that Cj is 
indeed a closed curve, and the definition of the permutation r shows that the right-hand side of (|36p 
does not depend on the choice of j G c €j. Equivalently, we can define Cj as the union over all j G 
of the curve segments {e^iX^) + iX'^9)) \ 9 G [0,vr]}. Let 

c A = du ...uc k c c , 

and let C' A C C be the bitangent set of Ca, namely the set of all line segments joining pairs of points 
of Ca at which Ca has the same tangent line. With these definitions, we have the following useful 
characterization of the singular set T,a- 

Proposition 5.1 [TTJ Theorem 3.5] S A = C A U C' A . 

Proof. According to Lemma \2.2\ Y>a is the set of all complex numbers of the form {Ax, x) where 
x G dW 1 is a normalized eigenvector of the Hermitian matrix H(9) for some 9 G [0, tt]. To describe 
that set, fix j G {1, . . . , m}, #o G R, and for 9 in a neighborhood of #o let 

z 3 {9) = (A Xj (9), Xj (9)) , 
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where Xj(9) is a normalized eigenvector of H{9) associated with the eigenvalue Xj(6) and depending 
smoothly on 9. General results in perturbation theory imply that such an eigenvector indeed exists 
[26] . Using the definition © of H(9), we find 

Re(e- i9 Zj (9)) = Re (e~ ie A Xj (9), Xj (9)) = (H{9)x j (9),x j (9)) = Xj(9) , 
Im(e-' w Zj (9)) = lm(e- id Ax j (9),x j (9)) = (H' (9)x j (9),x j (9)) = X'j(9) , 

since {H{9)x' j {9),x j {9)) + (H(9)x j (9), x' j (9)) = 2X j {9)Re(x j (9),x' j (9)) = due to the normalization 
condition. Thus 

Zj (9) = e ie (X j (9) + iX' j (9)) , hence z'^9) = ie i9 (Xj(9) + X" j {9)) . (37) 

These relations show that the curve Cj = {zj(9) \ 9 G M} is tangent, at each point Zj(9), to the straight 
line 

Lj(9) = {e ie (Xj(9) + ia) \ a G M} = {z G C | Re (ze- id ) = Xj(9)} . (38) 

In other words Cj is the envelope of the family of straight lines Lj(9), for any j G ^j. Since Cj C 
by construction, we have shown that T,a contains the curve Ca = C\ U . . . U Ck- 

However, it is important to realize that can be larger than Ca if the crossing set O C [0, 7r) 
defined above is nonempty. Indeed, assume that Xj(9o) = X p (9q) for some 9q G and some 
j,p G {l,...,m} with j 7^ p. For 9 in a neighborhood of #0) let Xj(9),x p (9) be smooth, normal- 
ized eigenvectors of H{9) corresponding to Xj(9),X p (9) respectively. Using the same notations as 
above, we have Zj(9o) ^ z p {9q) in general, because A^-(#o) 7^ ^p(#o)- Now, since (xj(9), x p {9)) = 
whenever Xj{9) ^ X p {9), we also have (xj(9o), x p (9q)) = by continuity. In particular, if a, (5 G C 
satisfy \a\ 2 + |/3| 2 = 1, then x = axj(9o) + /3x p (9q) is a normalized eigenvector of H(9q), and a direct 
calculation yields 

lm{e- Wo Ax,x) = {H'(9 )x,x) = \a\ 2 X'^q) + \/3\ 2 X' p (9 ) , 

whereas Re (e~ ie °Ax, x) = (H(9q)x,x) = Xj(9o) = X p (9q). This shows that Ha contains the line 
segment [zj(9o), z p (9o)], which by construction is tangent to the curve Cj at Zj(9o) and to the curve 
Cp at z p {9q). Repeating the same argument for all eigenvalue crossings, we conclude that Y^a contains 
the whole bitangent set C' A . Finally, it is clear from Lemma 12.21 that all points of Y<a either belong to 
Ca or to C' A . m 

Examples: 

1. ( The generic case) As is well known, within the space of all Hermitian matrices of size n > 2, the set 
of matrices having a multiple eigenvalue is a finite union of submanifolds of codimension at least three. 
This implies that, for a generic matrix A G M n (C), the Hermitian matrices H{9) defined by (|2|) will 
have simple eigenvalues for all 9 G 5 1 [IB] . In that situation, we denote by X\{9) < X%(9) < . . . < X n (9) 
the eigenvalues of H{9). Using the notations introduced above, we have m = n and the permutation 
r defined in ([35]) is simply 

A 2 ... n 
\n n — 1 ... 1 

If n = 2k is even, r has k cycles with i\ = . . . = 2; if n = 2/c — 1 is odd, r has one fixed point and 
k — 1 cycles of length 2. In all cases, the multiplicities mi, . . . , are all equal to 1, and the bitangent 
set C' A is empty by assumption. Thus £a = Ca is the union of k = [^4^-] closed curves. It is not 
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difficult to prove that the curve C\ associated with the cycle ^\ = (1 n) is smooth, strictly convex, 
and contains all the other curves C2, • • • , in its interior [19] . In particular, C\ = 91^(^4). Consider 
now the curve Cj = {zj(9) \ 8 £ R} for some j ^ l,n. If 5j(9) = Xj(0) + AJ(6>) is not identically zero, 
the formulas ([37]) show that the curvature of Cj at any regular point Zj(9) is strictly positive: 

Kj[6) = gM = |A#) + A?(6>)| > ° • 

This means that the tangent vector z'j(9) always rotates counterclockwise when 9 is increased. Never- 
theless, the whole curve Cj is not convex in general, because it may have a finite number of singular 
points corresponding to zeros of 6j(9). As is easily verified, simple zeros of 5j(9) correspond to cusp 
points of the curve Cj, see e.g. Fig. Q] where a generic example with n = 3 is represented. On the 
other hand, if 5j(9) vanishes identically, then z'J9) = and the curve Cj reduces to a single point. 
Under our generic assumptions, this can happen only if n is odd and j = (n + l)/2. As an example, 
the singular set of the matrix ^3 defined in (l33j) includes the isolated point {0}. 

2. ( The normal case) In contrast with the previous example, we now consider the particular case where 
the matrix A is normal. If Ai, . . . , X m denote the distinct eigenvalues of A, it is straightforward to verify 
that the eigenvalues of the Hermitian matrices H(9) defined in ([2]) are simply Xj(9) = Re(Aje ), 
j = 1, . . . , m. In view of (I37|) . this means that Zj(9) = Xj for all 9 G R, hence the curve Cj is reduced 
to the single point {Xj} for all j = 1, ... , m. Needless to say, the permutation r defined by (|35|) 
is the identity here. It follows that Ca = o~{A) = {Ai,...,A m }, and proceeding as in the proof of 
Proposition 15 . 1 1 we easily see that C' A is the set of all line segments joigning pairs of eigenvalues of A. 
This is in full agreement with the conclusions of Section 13.31 

Remark: We have seen that the curve Ca which generates the singular set Y<a is the envelope of the 
family of straight lines Lj{9) defined in (|38p . This geometric construction can be formulated in an 
equivalent way pQ, which is more conceptual and worth mentioning here. Assume for simplicity that 
the matrix A £ M ra (C) is "generic" in the sense of Example 1 above, and let a(£) be the homogeneous 
polynomial of degree n defined by 

a(0 = det(z I n + ZiAi + 6^2) , £ = (Co, 6, 6) G R 3 , 

where Ai,A2 are as in ([2]). Since a(£) has real coefficients, the equation a(£) = defines, in homoge- 
neous coordinates, an algebraic curve T in the projective plane RP 2 . Moreover, our genericity assump- 
tion on A implies that this curve is nonsingular: for each £ £ T, the tangent line to T at £ is uniquely 
defined and satisfies, in homogeneous coordinates, an equation of the form xo£o + + ^2^2 = 0. The 
set of all x = (xo, xi, X2) obtained in this way is again a curve V in RP 2 (the dual curve of T) given by 
the equation b(x) = for some homogeneous polynomial b. As is shown in [9], the degree of b does not 
exceed n{n — 1) if n > 2. Now, it is rather straightforward to verify that the curve Ca defined as the 
envelope of the family of straight lines (|38p is nothing but the restriction of the projective curve V to 
the subspace a^o = 1> namely z = x + iy 6 Ca if and only if 6(1, x, y) = 0. Thus Ca is a real algebraic 
curve in C ~ R 2 of degree at most n(n — 1) if n > 2. In the language of partial differential equations, 
the algebraic variety T is the characteristic variety of the symmetric hyperbolic system (|4|), and we 
shall see in Section [6] that the dual variety T' is related to the singular support of the fundamental 
solution of ([4|). Note that our genericity assumption on A precisely means that system (j4]) is strictly 
hyperbolic. 
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We conclude this section with a brief discussion of the number of tangent lines to the algebraic 
curve Ca which can be drawn from a given point. We recall that Ca = C\ U . . . U where each 
Cj is a closed curve associated with the cycle of the permutation ([35]) . For all J G {1, . . . , k} and 
all z G C \ Cj, we denote by Nj(z) the number of straight lines that are tangent to the curve Cj 
and contain the point z. Note that, since Cj was itself defined as the envelope of a family of straight 
lines, the tangent line to Cj is well defined even at singular points. In the degenerate case where 
Cj reduces to a single point {zj}, the set of tangent lines should be understood as the pencil of all 
straight lines through zj. Now, if z £ C \ Ca, we denote by N(z) the total number of tangents to the 
curve Ca = Ci U . . . U that can be drawn from the point z, with multiplicities taken into account: 

N(z) = m 1 N 1 (z) + ---+m k N k (z) , zEC\C A . (39) 

The following elementary properties of Nj(z) and N(z) will be useful. 

Proposition 5.2 For each J G {1, . . . , A;}, the number Nj{z) is constant in each connected component 
ofC\Cj. Moreover N(z) < n for all z G C\C A . 

Proof. Fix J G {1, . . . ,k} and pick j G ffj. For any z G C \ Cj, Nj(z) is the number of zeros of the 
function 

fj{9,z) = X^-Reize-' 9 ) 

for 9 in the interval [0,£jir). When z is varied, this number can only change if fj(9,z) has a double 
zero for some 6, but this would mean that z = Zj{9) = e ld (Xj{9) + iX'j(9)) G Cj, thus contradicting 
our assumption. Therefore Nj{z) is necessarily constant in each connected component of C \ Cj. On 
the other hand, we have the identity 

k 

det(iT(0)-Re(*e- i0 )In) = II II h^^T 3 ■ (^) 

J=ijeVj 

Fix z G C \ Ca, and consider both sides of (|40p as functions of 9. If multiplicities are taken into 
account, the number of zeros of the right-hand side on the interval [0,ir) is precisely N(z). But the 
left-hand side, being a trigonometric polynomial of degree at most n, cannot have more than n zeros 
in [0,7r). This proves the claim. ■ 

Remark 5.3 As was mentioned in the proof of Proposition 15, il if z G C \ Ca it is possible that 
fj(8,z) = f p (9,z) = for some 9 G [0,7r) and some j ^ p. This is the case, in particular, whenever 
z G C' A . However, if z E C \ T,a, we have 

d e fj{9, z) ■ d e f p (9, z) > whenever fj(9, z) = f p (9, z) = . 

Indeed, replacing A with A — zl n , we can assume without loss of generality that z = 0. // Xj(9) = 
X p (9) = for some 9 G [0, ir) and some j ^ p, then X'j(9) and X' p (9) have necessarily the same sign, 
otherwise (|37|) would imply that the origin belongs to the line segment [zj(9), z p (9)\, thus contradicting 
our assumption that z G" C' a . 

Proposition 15.21 asserts that the integer Nj{z) can only change if z crosses the curve Cj. In fact, 
if the crossing occurs at a regular point z G Cj, it is not difficult to verify that Nj(z) is decreased by 
two units if z crosses Cj from the convex to the concave side (i.e., in the direction of the local center 
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Figure 2: The singular set T,a (left) and the eigenvalues Ai (0), . . . , Ag(0) (right) are represented for a typical 
matrix A E M^CR). The singular set consists of three closed curves, one of which (the boundary of T4^(A)) is 
smooth, and the other two have cusps. In this example, the set Ha defined in (|46|) has six connected components. 

of curvature), and increased by two units if z crosses Cj in the converse direction, see [H Section 4.1] 
or Section 16.21 below. These simple rules give an efficient algorithmic way to compute Nj(z), and 
hence N(z), in concrete examples. Consider for instance Fig. [21 where the singular set T,a of a generic 
matrix A G Ms(lR) is represented. Here XU = Ca = C\ U C2 U C3, where C\ is the boundary of W(A), 
C2 is the closed curve with two swallowtails, and C3 is the triangle with three cusps. The set C \ 
has 11 connected components, on which N(z) is equal to 5, 3, or 1. Using the crossing rules above, it 
is easy to verify that N(z) = 5 in six different regions: inside both swallowtails, inside the three tips 
of the triangle, and outside W(-A). Such a result is definitely more cumbersome to obtain by counting 
directly the number of tangents to Ca from a given point. 

6 Qualitative properties of the numerical density 

Equipped with the results of Sections 0] and El we now derive some of the main properties of the 
numerical density Ja of an arbitrary matrix A G M n (C). We first establish an explicit formula for the 
derivatives of order n — 2, which allows us to prove that /a is polynomial in some distinguished regions 
of the complex plane which can be characterized geometrically. We next show that, for a generic matrix 
A G M n (C), the density Ja is of class C n ~ 3 if n > 3. Finally, as announced in the introduction, we 
prove that the fundamental solution of the linear hyperbolic system (jU) can be represented in terms 
of derivatives of the numerical density /a- In particular, the lacunas of system (j4]) are precisely the 
polynomial regions described in Corollary 16.21 

6.1 Polynomial regions 

Most of what we know about the numerical density of nonnormal matrices is based on the following 
result: 

Proposition 6.1 Let A G M„(C) ; where n > 2, and let V = P(d x ,d y ) be a homogeneous differential 
operator of degree n — 2. Then, for all z = x + iy G C \ Y,a, we have 




(41) 
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where 

A(9,z) = det( H(9) - Re{ze- ie )I n ) , 9 G S 1 . (42) 



In (|4ip the symbol f.p. denotes the finite part of the integral in the sense of Hadamard, but in 
many situations it is sufficient to take simply the Cauchy principal value, as in (|24|) . This is the 
case in particular for generic matrices in the sense of Section [5] (Example 1), because all multiplicities 
mi, . . . , rrifc are all equal to one and formula (|40p then shows that the map 9 i— > A(9, z) has only simple 
zeros on S 1 if z G C \ XU- 

In fact, using the analyticity of the integrand, it is possible to rewrite (|4ip in a slightly different 
form which is appropriate for further analysis. As in Section 14.11 we set w = e 2ie and we observe that 

P(cos6>,sin6>) = w x ~%P(w) , A(9,z) = w~^A(w,z) , 

where (if z = x + iy) 

PH = P(— , A(w,z) = det(- + — -(*— + y— )l„ 



Thus (|4Tj) is equivalent to 



4vH /| u |=iA(w,2) 

i(n-l)! / P(u>) , i(n-l)\ f P(w) , 

1 ; ■ dw + V o / * -/ \ du> , (43) 



8vr 2 7| tu | = i_ e A(w,z) 8vr 2 /j^i+g A(w, z) 

where e is any sufficiently small positive number, depending on z. In the particular case where P = 1 
and A G M2(C) is given by ([23]) . we recover (f25|h 

Proof of Proposition 16.11 We first consider the generic situation where the Hermitian matrices 
H(9) have simple eigenvalues Xi(6) < A2 (#)<••• < A n (#) for all 9 G 5 . In that case, the £>-spline 
representing the numerical density of ^h{9) can be given an explicit expression using the divided 
difference formula (fT4"P . (fT6|) : 



B[AiW,...,A n (e)]( S ) = (n-l)E Tr ^- ) S)r ' 



, ,. M-\ k {0)) 



In particular, differentiating (n — 1) times with respect to s, we find 

^■^■•■•■^iw = ("-^t nj(^w) - (44) 

Let = P(d x ,d y ) be a homogeneous differential operator of degree n — 2. From the representation 
formula ([22]) . we deduce at least formally 

(Vf A )(x + iy) = -*- / P(cos^,sin0)-H5 (n - 1) [Ai(e),...,A n (^)](xcos0 + ysin6')de . (45) 
47T J 5 i 
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To evaluate the integrand in the right-hand side, we start from (|44p and recall that the Hilbert 
transform (with respect to the variable s) of the Dirac measure S(s — A) is the distribution - p.v.-^r. 
We also use the identity 



it it 

y t = TT- 



which holds for any collection of pairwise distinct nonzero complex numbers /ii, . . . , fj, n . We thus find 

« b <->[, iW ,...,a„«,) 1W . _(^ P ,.n^_ . _(^, v ._i_ . 

Setting s = x cos # + y sin and inserting this expression into (|4"5"j) , we obtain (|4"Tj) . 

The calculations so far are formal, but they can be justified if we assume that z = x + iy £ C\ T*a- 
In that case, we know from the proof of Proposition 15.21 that the map 9 1— > A(6,z) has only simple 
zeros on S , because this is the case for each of the factors fj(9,z) in (|40p and, by assumption, the 
eigenvalues Xj(0) are all distinct for 9 £ S . It follows that the integral in (|4ip is well-defined in the 
sense of Cauchy's principal value, and depends smoothly on z G C \ S^. This in turn implies that the 
density Ja given by (f22l) is smooth on C \ T,a, as it should be, and that the calculations above are 
correct. 

To conclude the proof of Proposition 16.11 it remains to verify that (|4ip or (|43p holds for an 
arbitrary matrix A £ M n (C). To do that, we first observe that the singular set C C is an 
upper-semicontinuous function of A in the sense that 

SCHb^a) = sup dist(z,S J 4) > . 

Moreover the numerical density Ja, together with its derivatives, depends continuously on A in C\Sa- 
These rather classical facts can be established using, for instance, the representation formula ([8|) (we 
omit the details). On the other hand, it is not difficult to verify that the right-hand side of (|43l) 
depends continuously on A for each z £ C \ ~Ea- The crucial point here is that the polynomial map 
w 1— > A(w,z) keeps the same number of zeros on the unit circle (counted with multiplicities) if the 
matrix A is slightly varied; in particular, we can choose the same e > in (|43p for all matrices in a 
neighborhood of A. This property can be established using the factorization (|40p . Remark 15.31 and 
general results for perturbations of eigenvalues of Hermitian matrices, see [18]. Now, given A £ M n (C), 
there exists a sequence {Ag}^ e pj of generic matrices converging to A as £ — > 00. If z € C\ E^, we know 
that equation (|43p holds for Ag if I is sufficiently large, hence taking the limit t — > 00 and using the 
continuity properties mentioned above we obtain the desired equality. ■ 

As a consequence of Proposition 16.1} we now establish an important property of the numerical 
density in the regions of the complex plane where the number N(z) defined in (|39p takes its maximal 
value n. 

Corollary 6.2 Given A £ M„(C), let 

U A = {z£ C\S A N(z) =n} . (46) 

If n > 3, the numerical density /a is polynomial of degree at most n — 3 on each connected component 
of Ha- If n = 2, then /a = on Ua- 
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Proof. We shall prove that VJa vanishes identically on 11^ for any homogeneous differential operator 
of order n — 2. Indeed, in view of (|40p . (|42p . the assumption z E 11,4 implies that the map i— )■ A(0, 2;) 
has exactly n zeros (counting multiplicities) on [0, 7r). Equivalently, the polynomial A (to, z) has exactly 
n zeros on the unit circle {\w\ = 1}. But since A(w,z) has degree n, this polynomial has no zeros 
outside the unit circle, and using Cauchy's theorem we conclude that the first integral in the last 
member of (I43p vanishes. The second integral is also zero, because the numerator is a polynomial of 
degree at most n — 2, while the denominator has degree exactly n, hence the integrand P(w)/A(w, z) 
decays at least like \w\~ 2 as \w\ — > 00. Thus VJa = on n^, and the conclusion follows. ■ 

Remark 6.3 If z = x + iy E o~(A), then the polynomial A(w,z) has degree strictly less than n, and 
it follows from the above proof that z ^ 1T4. Thus cr(A) n 11^ = 0. 

The set 11^ is never empty, because it always contains the complement of the numerical range 
W(.A), where the density fA vanishes identically [261 H]. Moreover, in many situations, one or several 
components of IT^ are contained in W(A), in which case Corollary 16.21 gives nontrivial informations 
on the numerical density. For instance, if A E M n (C) is a normal matrix whose numerical range has 
nonempty interior, then N(z) = n for all z E C \ o~(A), and it follows from Corollary 16.21 that fA 
is polynomial of degree at most n — 3 in each connected component of C \ T,a, in agreement with 
Proposition 13.21 In the same spirit, if A = A% @ A2 with A\ E M ni (C) and A2 E M n2 (C), it is easily 
verified that N(z) = n = n\ + ri2 if z £ W(Ai) U W(A2), hence fA is piecewise polynomial outside 
W(^4i) U W(A2). Finally, Fig. [1] shows a typical example of a matrix A E M3(C) for which Ha has a 
component inside W(A), on which the density fA is identically constant by Corollary 16.21 

6.2 Generic regularity results 

Our purpose here is to establish regularity results for the numerical density fA in the whole complex 
plane, and not only outside the singular set T,a- For simplicity, we assume henceforth that our matrix 
A E M n (C) enjoys the following (generic) properties: 

HI: The eigenvalues of © satisfy Ai(0) < . . . < A n (0) for all 9 E S 1 . 

H2: For all j E {1, . . . , n}, the function Xj(9) + A"(0) is not identically zero. 

The second assumption guarantees that the curve (f3"6|) associated with Xj is not reduced to a 
single point. This of course is possible only if n > 2, and the two-dimensional case n = 2 is completely 
treated in Section l4.1i So we can assume that n > 3, and we have the following result: 

Proposition 6.4 Assume that A E M n (C) satisfies HI, H2 above. If n > 3, the numerical density 
f A : C is of class C n ~ 3 . 

We have already seen that both hypotheses HI, H2 are necessary, in general, for the conclusion 
of Proposition 16.41 to hold. For instance, if A E M^(C) is a normal matrix whose numerical range 
W(A) has nonempty interior, the numerical density fA is proportional to the characteristic function 
of W(A) and is therefore discontinuous on dW(A). A more subtle example is provided by the matrix 
A 3 defined in (33]): if \a\ 2 + \b\ 2 = 4, we have here Ai(0) = -1, A 2 (0) = 0, A 3 (0) = 1 for all E S 1 . 
Thus HI is satisfied, but obviously not H2, and the explicit formula (|34p shows that the numerical 
density of A3 is discontinuous at the origin. 
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Proof of Proposition E3J Assume that n > 3 and that A G M n (C) satisfies HI, H2. If Q = 
Q(d x ,d y ) is a homogeneous differential operator of degree n — 3, we have as in ([15]) : 

(Q/a)(z + %) = / Q(cos0,sin0)^s( n - 2 )[Ai(0),...,A„(^)](xcose + ysin^)d0 , (47) 

where 

Here H : R — > [0, 1] denotes the Heaviside function. Taking the Hilbert transform of (|48p with respect 
to s and using ([37]) . we thus find 

/■^ , w \ (n— 1)! /" Q(cos#,sin#) , , „ .„ . , , „ 

n . 

= V / ^-(5) log 1/^,^)1 dfl, (49) 

where fj(0, z) = Xj(9) — Re(ze~ l9 ) and hj : S 1 — > R is a smooth function. Note that the integrand in 
(f4~9~|) is 27r-periodic, because by assumption HI this is the case for all eigenvalues Xi(9), . . . , A n (0). 

It remains to show that each integral in the right-hand side of (|49p defines a continuous function 
of z £ C. Replacing A by A — zl n (an operation which does not affect the properties HI, H2), we see 
that it is sufficient to prove continuity at z = 0. This in turn is obvious if Xj(0) = fj(9,Q) does not 
vanish, so from now on we focus on the case where Xj(6) has (isolated) zeros on S . Using a partition 
of unity, we can treat each zero separately, so it is sufficient to consider the case where Xj(8) has a 
single zero of order q > 1 at 6 = 0, and hj{6) is localized near the origin. By analyticity, for z close 
to zero we have the factorization 

i 

fi&z) = g(9, z) H(9-fi p (z)) , 
P =i 

where the (possibly complex) roots /j, p (z) depend continuously on z, with /x p (0) = 0, and g(-,z) does 
not vanish in a neighborhood of zero. The quantity we have to study is therefore 

V / fc i (0)log|0-A» p (z)|d0+ / ^(6)^(9, z)\d9 . 

The last term is clearly a continuous function of z. In the integral involving fx p , we make the change 
of variables 9 = t + ~Re/i p (z) and observe that log |i| < log |i — ilmii p (z)\ < in a neighborhood 
of (t, z) = (0,0), hence continuity with respect to z follows from Lebesgue's dominated convergence 
theorem. ■ 

Once continuity of the derivatives of order n — 3 has been established, we can obtain further 
regularity results for the numerical density by using the representation formula (I4ip or (I43p . In view 
of Lemma 12.31 it is sufficient to study the density in a neighborhood of a point z £ S^. Since 
Y,a = Ca = C\ U . . . U Ck by assumption HI, there exists j S {1, . . . , n} such that z G Cj, and for 
simplicity we also assume that z £ C p for all p ^ j. This means that the function 9 t— > fj(9, z) has a 
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zero of order m > 2 at some point G 5 , and only simple zeros for 9^9 mod 7r; moreover, if p 7^ j, 
1 — ^ fp{0,z) has only simple zeros and does not vanish for 9 = 9. 

Without loss of generality, we assume from now on that z = and (9 = 0, and we first consider 
the simplest case where m = 2. This means that Xj{9) = ^9 2 + 0(9 S ) near 9 = 0, for some If 
z = x+iy is sufficiently small, it follows that the analytic function 9 1— > fj(9) = Xj(6) — (x cos 9 +y sin 9) 
has exactly two zeros 9±(z) in a neighborhood of the origin, which satisfy 

0±{z) = ^(y±sign(a)vV + 2as) +0{\x\ + \y\ 2 ) . 

The critical curve Cj, which is the set of all points z for which 0+(z) = 9-{z), is therefore given by 
the equation x = —-^y 2 + C(y 3 ) m a small ball B around the origin. Moreover B \ Cj = B r U B c , 
where B r is the set of all z G B for which the roots 9±{z) are real and distinct, whereas z G B c when 
9±(z) are complex conjugate with nonzero imaginary part. As is easily verified, the local center of 
curvature of Cj is located on the side of B c . 

Now, let V = P(d x , d y ) be a homogeneous differential operator of degree n — 2, and consider the 
expression of (VfA)(z) given by (|43p . Assume that, in the right-hand side, the parameter e > is 
chosen in such a way that the slit annulus A t = {w G C|0 < |1 — \w\\ < e} contains the points 
w±(z) = e 2l9 ±( z ) for all z € B c , but that Aie does not contain any other root of the determinant 
A(w,z) for z G B (these conditions are easily achieved by choosing first e and then B sufficiently 
small). Under these assumptions, the right-hand side of (|43p defines a smooth function of z G B, 
which coincides with {T > fA){z) if z G B r but not if z G B c . Indeed, in the latter case, we have 
to consider in addition the roots w±(z) of A(w,z) which are not taken into account by the fixed 
integration contours in (I43j) since |1 — |iu±(.z)|| < e. Using Cauchy's theorem, we easily obtain 

(Vf A )(z) = (Vf A ) ies (z) + — w _ {z) _ w+{z) • 

where (T'fA)reg(z) denotes the regular part of (VfA)(z), given by the right-hand side of (|4"3"j) with 
fixed e, and 

P(w)(w — w + (z))(w — W-(z)) 



Q(w, z 

In particular, we have for all z G B c 
(Vf A )(z)-(Vf A ) ieg (z) 



A(w,z) 



< C - < C - < ° 



\w+{z)-w-(z)\ - \0+{z) - 0-{z)\ ~ dist^Cj) 1 / 2 



and this estimate is sharp if P(l) = P(1,0) / 0, because in that case Q(w,z) does not vanish near 
w = 1 if z G B. Summarizing, we have reached the following important conclusion: Near a regular 
point z G Cj of the critical set T, A of a generic matrix A G M n (C), the derivatives of order n — 2 of the 
numerical density f A are smooth on the convex side of the curve Cj, and blow up like dist(z, Cj)~ l l 2 
on the concave side. If n = 2, this is in full agreement with the explicit formula (|27p obtained in 
Section [4.1i If n > 3, we deduce after integrating that the numerical density f A is of class C n_5 / 2 in 
a neighborhood of such a point z. 

Using the same techniques, it is also possible to study the regularity of the numerical density near 
more singular points z G T, A - On the typical example represented in Fig. [21 we see that the following 
two cases have to be analyzed: 
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i) Crossings, which occur when two critical curves Cj and C p intersect transversally at z. This situation 
can be treated exactly as before, except that one has to consider four distinct regions near z, instead 
of two. 

ii) Cusps, which arise whenever one of the functions \j{9) + A''(0) has a simple zero. Here we can 
repeat the analysis above, assuming that Xj(0) = ^9^ + 0{9 4 ) near 9 = 0. In a neighborhood of 
z = 0, we find that fj(9, z) has either three real roots (for z inside a cuspidal domain with tip at z), 
or one real and two complex conjugate roots (outside the cusp). Using the same argument as before, 
we conclude that V/a is smooth inside the cusp, but blows up on the other side of the critical curve 
Cj. Altogether, the numerical density /a is of class (7 n ~ 8 / 3 near the cusp, see [26j Section 4.3] for a 
similar analysis of the singularities of the fundamental solution of (|4|). 

Under generic assumptions on the matrix A, all intersections are transversal and all cusps are 
non degenerate, so that the singular set Ca is a generic curve in the sense of real algebraic geometry, 
and the singularities of the numerical density /a can be completely analyzed using the techniques 
described bove. In particular, the derivatives of order n — 2 of the numerical density Ja are locally 
integrable, and since we know from Proposition 16.41 that Ja £ C n_3 we conclude that the relation (|4ip 
holds everywhere (in the sense of distributions), and not only on the complement of S^. 

6.3 Connexion with the fundamental solution 

As was explained in the introduction, the numerical measure of a matrix A 6 M n (C) is related to the 
fundamental solution of the hyperbolic system This connexion can be established rigorously by 
comparing the expression (|4ip for the derivatives of the numerical density /a with the representation 
formulas for the fundamental solution E{t,x) of (j3|), which can be found e.g. in [4"1 126j. 

In what follows, we identify C with M 2 , and we denote by x = (x\, X2) the points of the Euclidean 
plane. With a slight abuse of notation, we write Ja(x) instead of Ja{z), and we consider as subsets 
of R 2 the various regions associated with A, such as IU(^4) or S^. To derive a representation formula 
for the fundamental solution E(t,x), we take the Radon transform of © with respect to 1 6R 2 , and 
obtain the one-dimensional hyperbolic system 

d t E(t,s,9) + H(9)d s E(t,s,9) = I n S t=0 ® 5 S=0 , 

where H(9) is given by ([2]) and E{t,s,9) denotes the Radon transform of E(t,x). Using the method 
of characteristics, we easily find 

11 

E(t,s,9) = ^25{s -t\ j {9))P j {9) , t>0, 
3=1 

where \\{9) , . . . , \ ri (9) are the eigenvalues of H{9) and P\{9) , . . . , P n (9) the corresponding spectral 
projections. Now, if we invert the Radon transform as in ([22]) and use the identity 

n 

we arrive at the representation formula 

E(t,x) = -^f.p.y {[x lC os9 + x 2 sin9)I n -tH{9)^ 2 d9 , t>0, (50) 
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which coincides with Eq. (4.4a) in [3]. Arguing as in Section [6.1( one can show that equation (|50p is 
rigorously satisfied for all (t, x) G (0, oo) x R 2 with | G R 2 \ T, A , and that E(t, x) is smooth in that 
region of space-time. 

To compare the numerical density with the fundamental solution, it is natural to extend f A to a 
homogeneous function of space and time by setting 

T A (t,x) = t n ~ 3 fA^) , t>0, xGM 2 . (51) 

We then have the following result: 

Proposition 6.5 For any A G M„(C), there exists a matrix-valued homogeneous polynomial Q of 
degree n — 1 such that 

E{t,x) = Q(d t ,d Xl ,d X2 )T A (t,x) , (52) 
for all (t, x) G (0, oo) x R 2 with f G R 2 \ E A . 

In particular, if Q is a connected component of the region n A C R 2 defined by (@SJ), we know from 
Corollary 16.21 that f A (x) is polynomial of degree at most n — 3 in fi, and (I5ip then shows that J- A (t, x) 
is also polynomial of degree at most n — 3 in the half-cone C+(0) = {(t, x) G (0, oo) x R 2 | | G f2}. By 
Proposition 16. 5[ we conclude that E(t,x) = in C + (£l). We thus have 

Corollary 6.6 The fundamental solution E(t,x) of a matrix A G M n (C) vanishes for all (t,x) G 
(0,oo) x R 2 such that | belongs to the region H A C R 2 defined by (JM|). 

In the language of partial differential equations, the domain of influence of the origin for system 
(gD is the half-cone C + (D) = {(t,x) G (0, oo) x R 2 | f G D}, where D C R 2 is the complement of 
the largest connected open region on which E* = E(l, •) vanishes. If in addition I?* = in some 
open region L C D, we say that L is a lacuna of the hyperbolic system (j3j). With this terminology, 
Corollary 16.61 asserts that each connected component of 11^ either lies outside the domain of influence 
of the origin, or is a lacuna. This important geometric characterization of lacunas is originally due 
to Petrovsky [22], and was thoroughly discussed in [U[26] and from a more algebraic point of view in 
[21 [3]. Strictly speaking, it gives only a sufficient condition for the occurence of lacunas, but further 
work allows to show that all stable lacunas satisfy this criterion. In other words, in generic situations, 
the converse of Corollary 16.61 also holds: each open region on which E* vanishes (in particular, any 
lacuna) belongs to IT4. 

As a typical example, consider the matrix A G M^(C) defined by (|55j) . whose numerical density 
is represented in Fig.[TJ The polynomial region II^ has just two components here: the exterior of 
D = W(A), where both f A and E* vanish, and the interior of the cuspidal triangle, where E* = and 
f A is identically constant, which is therefore a lacuna. A more complicated situation is depicted in 
Fig. [21 where the polynomial region 1T4 has six connected components, among which five correspond 
to lacunas. In a different spirit, it is also interesting to consider the nongeneric matrix (I57p whose 
numerical density is studied in Section [71 The domain of influence of the origin is the union of the 
closed unit disk and a single point {a}, so the numerical range Vl^(A) = conv(L>) is substantially larger 
than D if \a\ > 1, see Fig.[5j In that case, the region 11^ has again two components, none of which is a 
lacuna: the exterior of W (A), and the interior of the triangular region W{A)\D where f A is identically 
constant. We see on these examples that the numerical density allows to distinguish between various 
regions where the fundamental solution vanishes identically, and which are nevertheless of rather 
different nature. 
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Before proving Proposition 16.5^ we briefly verify its conclusion on a simple example. If A £ M 2 (C) 
is defined by (|3"T|) . system ^ reduces to 

dfU\ + d Xl u 2 - id X2 u 2 = , d t u 2 + d Xl ui + id X2 u\ = . 

Combining both equations, one verifies that dfv,j = Auj for j = 1, 2, hence the fundamental solution 
E(t, x) can easily be computed using Poisson's formula for the solution of the wave equation in two 
dimensions. In agreement with Proposition 16.51 the result is: 

where according to ([32]) . (f5"T|) 



Proof of Proposition 16.51 If n = 1 the conclusion is trivial, because both E* and Ja vanish 
identically outside Y<a (which is reduced to a single point), so we assume henceforth that n > 2. 
If "P = P(dt,d Xl ,d X2 ) is a homogeneous differential operator of degree n — 2, then using (|51~j) it is 
straightforward to verify that 



P 



(53) 



whenever j £ M?\T,a- Here £ S M 2 denotes the argument of the function /a, which has to be replaced 



by j after differentiation. We warn the reader that equality 
Applying Proposition 16. 1\ we deduce that 



holds only if P is of degree n — 2. 



P{d t ,d Xl ,d X2 )FA{t,x) 



(n-1)! 
4vr 2 t 



f.p. 



P(— f • e<9,cos#,sin#) 
det(P(#) - f 



in 



d0 



(54) 



where eg = (cos 0, sin 9). On the other hand, starting from f)50f) . we observe that 



d_ 

ds 



x 



s + ~- eg)I n - H{9) 



-i 



s=0 



By Cramer's rule, the inverse of the matrix SI n — H{9) = SI n — A\ cos 9 — A 2 sin 9, with S = s + 
has the following form 



SI n ~ H{9) 



1 



A(0,S) 



S Pq(S, cos#,sin#) + cos9 Pi(S, cos#,sin#) + sin#P2(5', cos#,sin/ 



where A(#, S 1 ) = det(H (9) — S I n ) and Po, -fi) ^2 are matrix-valued homogeneous polynomials of degree 
n — 2. The idea is now to insert this expansion into the right-hand side of ([50]) and to use ([54")) to 
express the result as a derivative of order n — 1 of the numerical density Fa- 
We begin with the term involving Pi, and remark that 



cos 9 



d (P\{s + f • eg, cos 6, sin0) 



9s 



A(0,8 + f 



9 f P\{j " e 0> c °s #j sin( 



s=0 



9xi 



A(0,f -e*) 
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The corresponding contribution to (|50p is thus 



1 9 r_ f Pi(j ■ e 9 ,cos9,sm6) jn 1 d 



Similarly, the term involving P 2 gives the contribution 

1 d 

E 2 (t,x) = - 7— jTj -Q^P2{-dt,d Xl ,d X2 )T A (t,x) . 

Finally, to treat the expression containing Pq, we observe that 

d f (s + f • ee)Po(s + § • eg, cos 0, sin 9)\ f P (j ■ ee,cos9,sm9) 



x 



ds V A(#,s + §-e e ) / s=o v ^\ A(9, 

and using (f53]h (f54"|) we obtain the contribution 

1 1 a 

£Jo(t,x) = -j^—^ t {l + x-V x )Po{-dt,d Xl ,d X2 )F A {t,x) = 7— jT| —Po{-d t ,d Xl ,d X2 )F A {t,x) . 

Recalling that E(t, x) = E (t, x) + Ei(t, x) + E 2 (t, x), we arrive at ([52]) . ■ 

7 Three-dimensional examples 

To illustrate the results of the previous sections, we consider here four concrete examples which, 
according to [19} Section 1.7], give a complete picture of what can happen for three-dimensional 
matrices. The two-dimensional case, which is much simpler, was already treated in Section 14.11 and 
the references [U HJ [16] include a detailed study of the singular set Yi A for a few higher-dimensional 
examples. 

Example 1. We first consider the 3x3 matrix 

A = ( -1 1 i J . (55, 




which is generic in the sense that hypotheses HI, H2 in Section T6.2I are fulfilled. The eigenvalues of the 
associated Hermitian matrix H(9) satisfy Xi(9) < X 2 (9) < ^3(9) for all 9 £ [0,n], as can be seen from 
Fig. [3] (right). The critical set T, A = C A = C\ U C 2 is an algebraic curve of degree 6 consisting of a 
smooth ovate curve C\ enclosing a cuspidal triangle C 2 (Fig.El left). The component C\ corresponds to 
the eigenvalues Xi(9), Xs(9) while C 2 is associated with X 2 {9}. All multiplicities are equal to one, and 
the number N(z) defined by fj39|) is equal to 3 outside C\ and inside C 2 , and to 1 in the intermediate 
region. The numerical density f A is continuous, identically constant inside C 2 , and vanishes outside 
C\. Moreover f A is Holder continuous with exponent 1/2 across C\ and C 2 , except at the cusps. The 
level lines of f A are represented in Fig.[TJ 

Example 2. We next consider a nongeneric matrix 

1\ 

(56) 
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Figure 3: The critical set £a (left) and the eigenvalues Xj(0) (right) are represented for the matrix ([55]). 



for which the critical set Y<a can be computed exactly. Indeed, if H{9) is the Hermitian matrix ([2]), it 
is easy to verify that det(A/3 — H(9)) = A 3 — |A — j cos 9, hence 

A lW = cos(^), X 2 (9) = cos(^) , A 3 (0) = cos(^) . 

Thus the permutation ([35]) is just a cycle r = (1 2 3), and applying ([37]) we easily find that the critical 
curve C A is the cardioid defined by C A = {\{2e i(t> + e 2 ^) \ <j> G S 1 }. Since Ai(0) = A 2 (0), the bitangent 
set C' A is not empty and consists of the line segment joining the points —1/2 ± i/(2\/3), see Fig.0] 
(left). Altogether we have £a = Ca U C' a , and we observe that Syi encloses a convex region of the 
complex plane which is of course the numerical range W(^4). The index N(z) defined by (139p is equal 
to 3 outside Ca and to 1 inside. The numerical density /a vanishes outside W(A) and is equal to a 
nonzero constant inside the cuspidal region, in agreement with Propositions 12.6 1 and 16. 11 In particular, 
/a is discontinuous along the line segment C' A . 



i.o- 




-1.0- 



Figure 4: The critical set Ha (left) and the eigenvalues \j(6) (right) are represented for the matrix ([56l) . 

Example 3. The matrices considered so far were unitarily irreducible. In contrast, the matrix 

/0 2 0\ 

A = , where a G C , (57) 
\0 a J 
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is the direct orthogonal sum of the two-dimensional Jordan block (|3ip and the one-dimensional matrix 
(a). The numerical density /a can therefore be computed using Proposition I2.7[ and without loss of 
generality we can assume that a > 0. However, we have to distinguish between three cases: 

i) If < a < 1, the numerical range VF(^4) is the closed unit disk, and the numerical density has the 
following expression: 



/(*) 



irVT 



: argch 



1 — az\ 



v/(l-a*i) 2 -(l-H 2 )(l 



\z\ < 1 , 



which reduces to (|34|) when a = 0. In particular /a vanishes on the unit circle, has a logarithmic 
singularity at the point {a}, and is otherwise smooth. The singular set T,a is the union of the unit 
circle and the point {a}. 

ii) In the limiting case a = 1, the numerical range is still the closed unit disk, but the formula 



/(*) 



iyr 

7T 



\z\ < 1 , 



1-21 ' 

shows that the numerical density has now an algebraic singularity at the boundary point z = 1. 

iii) When a > 1, the numerical range VF(A) is the convex hull of the union of the unit disk and the 
exterior point {a}. Within this region, the numerical density satisfies 



/(*) 



1 



arccos 



az\ 



+ (1-M 2 )(a 2 



when \z\ < 1 



and f(z) = 1/V a 2 — 1 when \z\ > 1. As is easily verified, the eigenvalues of the Hermitian matrix 
H(9) are Xi(6) = —1, A2(6>) = 1, and \3(0) = acosO, see Fig. [5] (right). The algebraic curve Ca 
consists of the unit circle (associated with Ai,A2) and the point {a} (corresponding to A3), but the 
bitangent set C' A is not empty and consists of two line segments, see Fig.0 (left). 





Figure 5: The critical set (left) and the eigenvalues \j(9) (right) are represented for the matrix (|57l) with 
a = 2. 



Example 4. As a final example, we consider the case of a normal matrix A G M^(C) whose eigenvalues 
Ai, A2, A3 are not colinear. Then the numerical range is the triangle with vertices {A1, A2, A3} and the 
numerical density is a multiple of the characteristic function of W(A). In that situation = C^UC^, 
where Ca is the set of all vertices and C' A the set of all edges of the triangle. 
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8 Statistical properties of the numerical measure 



In this section, we study the numerical measure from a statistical point of view, and we establish various 
convergence results which show that the measure [i A of a large matrix A 6 M n (C) is concentrated in 
a neighborhood of the barycenter of the spectrum cr(A). 



8.1 Concentration phenomena 

Proposition 8.1 For any A £ M n (C), the first moment of the probability measure fi A is the normal- 
ized trace of A: 

Ji A = zdfi A (z) = - Tr(A) , 
Jc n 

and the variance of \x A is given by 

Vax(p A ) = f \z-Jl A \ 2 df, A (z) = — !— (- Tr(A*A) - - Tr(A) *) . (58) 

Proof. Let a,& denote the entries of A. Applying definition ([1]) with <j)(z) = z, we have to compute 
the average of YLjk a jk x ~j x k over the unit sphere. By symmetry, the average of x^x^ is equal to zero if 
j 7^ k and to 1/n if j = k. Thus 

~Pa = -( a n H ha„ n ) = - Tr(yl) . 

n n 

To compute the second moment, we take <j>(z) = \z\ 2 and proceed in exactly the same way. By 
symmetry, the average of Tjx^xlx m is zero unless j = k and I = m, or j = m and I = k. Moreover, it 
is easy to verify that the average of |xj| 2 |x£| 2 is equal to r n if j = £ and to s n if j ^ £, where 



/ |xi| 4 dcr(2;) = 2 , s n = I \x 1 \ 2 \x 2 \ 2 da(x) 
Jem™ n{n+l) J m n 



1 



nyn + 1) 



Thus 



\z\ 2 &n A {z) 

c 



j,k,£,m 
n 



= Sn)X a i3 a M + Wj?) + r n \ a jj\ 2 
n 

This gives the desired result, since Var(^) = J c \z\ 2 dfi A (z) — \ji A \ 2 . ■ 

Now we consider a sequence of matrices {A n } n >i such that A n £ M n (C) for each n > 1. As is 
well known, we have 

- Ti(A* n A n ) < \\A n \\ 2 < Tr«,4 n ) , 
n 

where \\A n \\ = sviY>{\\A n x\\ \ x G dW 1 }. As a consequence, if we suppose that ||^4 n || 2 = o(n) as n — > oo, 
it follows from (|58p that the variance of /U^ converges to zero as n — > oo. This gives: 
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Corollary 8.2 Assume that A n G M„(C) and that \\A n \\ 2 /n — > as n — > oo. Then the measure 
HA n — 5-p_ A converges weakly to zero as n —> oo. 

We recall that the numerical radius of a matrix A G M n (C) is defined by 

R(A) = sup{\z\ | z G W{A)} = sup{| {Ax, x)\\x G dB n } , 

and satisfies R(A) < \\A\\ < 2R(A) [15J. Thus, a sequence of matrices A n G M n (C) is uniformly 
bounded (in the operator norm) if and only if the numerical ranges VF(^4 n ) are all contained in a 
bounded region of the complex plane. Under this assumption, Proposition 18. II shows that the variance 
of HA n is 0{l/n) as n — > oo, so that the numerical measure is asymptotically concentrated in a disk 
of radius 0(1/ y/n) around the mean ~p-A n - 

In the introduction, we have observed that the numerical measure (j,a„ is the distribution of the 
random variable (A n X n , X n ) G C when the vector X n is uniformly distributed on the unit sphere dM n . 
With this interpretation, Corollary 18.21 is reminiscent of the weak law of large numbers in probability 
theory. Under slightly stronger assumptions, it is also possible to obtain a pointwise convergence result 
in the spirit of the strong law of large numbers. Without loss of generality, we assume from now on 
that Ti(A n ) = for all n > 1, so that the measure \iA n is centered at the origin. 

Proposition 8.3 Assume that A n G M n (C) satisfies Tr(^4 n ) = for all n > 1 and 

(logn) \\A n \\ 

sup ye < 00 • ( 59 ) 

n>l fl ' 

If for each n > 1 the random variable X n is uniformy distributed on the unit sphere dW 1 , then 
{A n X n , X n ) converges almost surely to zero as n — )• oo. 

Proof. It is clearly sufficient to prove the result for Hermitian matrices A n , because the general 
case then follows by considering the real and imaginary parts of {A n X n , X n ) . We thus assume that 
A n = A* n for all n > 1, and we denote by A n> i, . . . ,X n . n the eigenvalues of A n . For each n > 1, 
let Y n> i, . . . , Y n ^ n be independent and identically distributed complex random variables with density 
function fy[z) = 7r _1 e _ ' z ' 2 , z G C. In particular, we have E[\Y n ^ m \ 2k ) = k\ for each k G N. Since the 
Euclidean measure on dW 1 is the projection on the unit sphere of the standard Gaussian measure in 
C n , we obtain a uniformly distributed random variable on <9B n by setting X n = U n Y n /\\Y n \\, where 
Y n = (Y n> i, ... , Y ntn ) T and U n G U n (C) is a unitary matrix such that U*A n U n = diag(A ni i, . . . , A ni „). 
Thus 

/AY V \ An,l|-^n,l| ~\~ ' ' ' ~\~ \n,,n \ Yn,n \ _ Pn 

\A n Ji n ,Ji n ) - 2 |2 - jr- , 

1 n 2 1 n 2 1 >" 

Rn — / , A n m \Y n m \ — \ A n m ( y^ m 1) , Q n — 



where 



Yn 



By the strong law of large numbers, the denominator Q n converges almost surely to 1 as n — > oo, 
hence it remains to show that the numerator P n converges almost surely to zero. But this follows from 
classical theorems on the limiting behavior of weighted sums of independent random variables, see 
[23]. Since for each n > 1 the random variables X n ^ m = \Y n ^ m \ 2 — 1 (1 < m < n) are independent, 
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have zero mean and finite second order moment, and since by (|59[) the coefficients a n ^ m = n 1 X n ,m 
satisfy 

C 

max |<z nr n,| < — 77^- , 

l<m<n n i ' z iogn 

the results of [23J Section 3] imply that P n = Ylm=i a n,m,X n ,m converges almost surely to zero as 
n —7- 00. ■ 

8.2 Central limit theorems 

The results established so far show that for a sequence of traceless matrices A n £ M n (C) the numerical 
measure ha„ tends to concentrate on the origin as n — > 00. Under stronger assumptions, we now prove 
that the rescaled measure /J*^A n converges to a Gaussian distribution, as in the classical central limit 
theorem. We first consider the Hermitian case, which is somewhat simpler. 

Proposition 8.4 Let {^4 n } n >i be a sequence of Hermitian matrices such that A n G M„(C) and 
Tr(A„) = 0. We assume that 

(logn)|L4 n || ^ Q ^ 1 ^ ^ cr^ > . (60) 

Then the rescaled numerical measure ^^A n converges weakly to the normal distribution N (0, a 2 ) as 
n — > 00. 

Proof. We use the same notations as in the proof of Proposition 18.31 For each n > 1, the numerical 
measure \xa u is the distribution of the random variable 

ry |^n,l I ~\~ ' ' ' ~\~ X n ^ n \Y nn \ P n 



l-in.ll t 1 \ 1 n,n\ ^cn 



where \ n ,l, . . . , X n ,n denote the eigenvalues of A n and Y n A, • • • , Y n<n are independent complex random 
variables with density function /y(z) = 7r _1 e - '^ . Since Q n converges almost surely to 1 as n — V 00, 
we have to show that n}l 2 P n converges in law to A/"(0, a 2 ). 

To do that, we use the Lindeb erg- Feller theorem for triangular arrays of random variables [8l 
Section 2.4.b]. Let 

X n ,m — 1= X>n,m (|^n,m| -Q > 1 — ^ — ^ i 



so that r^l 2 P n — X n \ + •••-)- X n ^ n . For each fixed n ^ 1, the random variables X nm are independent 
and satisfy E(X n ^ m ) = for m = 1, . . . , n. Moreover, 



£ £(|X n ,J 2 ) = - £ A 2 , m = - Tr(^) — a 2 > . 

m=l m=l 

Finally, for any e > 0, we have 

E(\X n , m \ 2 ; |X n , m | = (\z\ 2 — l) 2 e-l-l 2 dz , 

* y J-Sn,.?~n..e 
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where -D n ,m,e = {z £ C | A^ m (|z| 2 — l) 2 > ne 2 }. Thus, using the first assumption in (f60]h we obtain by 
a direct calculation 

n 

y~] E(\X nm \ 2 ; \X nm \ > e) < n sup E(\X nm \ 2 ; \X nm \ > e) > . 

~ 1 l<m<n n ~>°° 

Invoking the Lindeberg-Feller theorem, we conclude that n 1//2 P n = X n \ + • • • + X n ^ n converges in law 
to Af(0, a 2 ), which is the desired result. ■ 

Remark. Since the numerical density of a Hermitian matrix is a l?-spline, Proposition 18.41 shows 
under very general assumptions that i?-splines of degree n satisfy a central limit theorem in the limit 
n — > oo. In the particular case of uniform i?-splines, this result was obtained by Unser et. aJ. in [25 1 . 

Before considering more general matrices, we would like to mention an alternative proof of Propo- 
sition EH which has its own interest. The starting point is a very nice formula for the moments of the 
numerical measure of a Hermitian matrix A £ M„(C). Fix k £ N and let Ai, . . . , A n £ M denote the 
eigenvalues of A. Using (fT71) or (fT8j) with (f>(x) = x k , we find 



/ x k dfJL A (x) = (n - 1)! / (tiAi + • • • + t n X n ) k dti 



dt n . 



where is the (n— l)-dimensional simplex defined in (I12p and t n = 1 — {t\-\ +t n -\). To evaluate 

the right-hand side, we apply the multinomial formula 

(x l + ... + x n ) k = V -x a , 

\a\=k 

where the sum runs over all multi-indices a £ N n of order \a\ = a% + ■ ■ • + a n = k. Here we use the 
standard notations X a = X" 1 ■ ■ ■ X" n and a! = (ai!) • • • (a«0- Now, it is not difficult to verify that 

a! 

t a dh... dtn-! = — , 

D„_! {n + k-iy. 

for any a £ N n with |a| = k. We thus obtain the following identity 

V ; \a\=k 

which shows that the k-th moment of the numerical measure \xa is the complete symmetric homoge- 
neous polynomial of degree k in the variables Ai, . . . , A n , divided by the combinatorial factor ( n+ t -1 ) 
which is just the number of terms in the sum. 

Using the Newton identities, the right-hand side of (|6ip can be decomposed into as a sum of 
products of elementary symmetric polynomials of the form pi = \\ + • • • + A^, see |2H Eq. (2.14')]. If 
we assume that Tr(^4) = 0, then p\ = and the number of nonzero terms in the sum is considerably 
reduced. Using these remarks, it is not difficult to show that, under the assumptions of Proposition[ 
the k-th. moment of the rescaled numerical measure ^^a„ satisfies 



if k is odd , 

2~k/2 a k_^d__ j£ £ ig even _ 
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Since the moments in the right-hand side are those of the normal law, we conclude that ^^A n 
converges weakly to M(0, a 2 ) as n — > oo. 

We now consider general matrices A n G M„(C), and obtain a central limit theorem by applying 
Proposition 18.41 to the Radon transform of A n . 



Theorem 8.5 Let {vl n }n>i be a sequence of matrices satisfying A n G M ra (C), Tr(A n ) = 0, and 
(log 7T.) ||^4tj. II /tt, 1 / 2 — > as n — >■ oo. We assume that 



1 



Tr(A* n A n ) 



-)■ a > , 



Tr(A 2 r 



-> b G [0, a) . 



(62) 



T/ien the rescaled numerical measure [i^A n converges weakly to the Gaussian measure f OQ (z)dz as 
n — > oo, where 



foo(x + iy) 



Tr\/a 2 — b 2 



g a-\-b a — b 



(63) 



Proof. For each 9 G S 1 , the Hermitian matrices H n (9) = i(e lt> A n + e l6 A^) satisfy the assumptions 



of Proposition 18.41 with 
a 2 = a{9) 2 



lim — Tr(A* n A n ) + Re(Tr(^)e 

n— »oo Zn 



-2i9- 



-(a + bcos(29)) . 



Let us denote by [x n and \x n $ the numerical measures of ^fri A n and ^ynH n (9), respectively. Since fx nt g 
is the Radon transform of /i n , the two-dimensional Fourier transform ft n (C) f° r £ = re% ® is precisely 
the one-dimensional Fourier transform of [i n fi evaluated at r |12j. Thus, applying Proposition 18.41 we 
find 

An(re ie ) = / e- ixr d^ e (x) ' 



e' lxr e 



-xy(2a(ef) dx 



-u(8) 2 r 2 /2 



o-(9)V2ir 

for any r > 0, 9 G S 1 . This shows that /i„ converges weakly as n — > oo to the measure on C 
defined by 

Inverting the Fourier transform, this gives d/ioo = foo(z) dz with as in ([63]) . ■ 



Remarks. 

1. If we assume for simplicity that ||^4 n || is uniformly bounded, we can suppose (up to extracting a 
subsequence) that ±-T?(A* n A n ) converges as n — > oo to some a > 0. However, we have to assume in 
(|62p that a > in order to get a universal Gaussian limit. 

2. Similarly, the first assumption in (|62|) implies that ^ Tr(^) converges, after extracting a subse- 
quence, to some b G C with |6| < a. Multiplying A n by a unit complex number, we can assume that 

< b < a, but in the borderline case where b = a the limiting measure foo(z) dz should be replaced 

(27ra)- 1 / 2 e - ;c2 /(2a) dx g, §y=0 . 



Example. Let us consider the Jordan block of size n: 

/0 1 



o\ 



1 
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Since Tr(A* A n ) = n — 1 and Tr(^) = 0, the assumptions of Theorem 18.51 are satisfied with a = 1 and 
6 = 0. Thus the numerical measure of ^JnA n converges weakly to the normal distribution ir~ 1 e~\ z \ dz 
as n — > oo. In this example, the measure HA n has in fact a radially symmetric density for all n > 1, 
see Section 14.21 

9 Perspectives 

As a conclusion, we briefly mention a natural extension of our work, which is left for future investi- 
gation. Recall that a homogeneous polynomial P G M.\Xq,X\, . . . , X d ] of total degree n is hyperbolic 
in the direction cq = (1, 0, . . . , 0) if, on the one hand, it has partial degree n with respect to the first 
variable Xq, and on the other hand the n roots of the univariate polynomial t i— > P(t,y) are real for 
every vector ?/ G R rf . Hyperbolic polynomials arise as principal symbols of hyperbolic differential op- 
erators of order n in d space variables, see As an example, if A\, . . . , A d G M n (C) are Hermitian 
matrices, the polynomial 

P{X Qi X x , ...,X d ) = det(A /„ - X 1 A 1 X d A d ) (64) 

is hyperbolic. In the particular case where d = 2, it has been conjectured in [20J, and proved in |14| . 
that all monic hyperbolic polynomials are of the form (I64p . This is no longer true if d > 3. 

It might be argued that a large part of our work is not really about matrices, but rather concerns 
hyperbolic polynomials in 2 + 1 variables. Indeed, if A G M n (C) and A\, A2 are as in ([2]), the 
eigenvalues \i(9), ■ ■ ■ , A n (0) of the Hermitian matrix H{6) are the solutions of the equation 

P4A,cos0,sin0) = , 6 G S 1 , 

where Pa(Xq, X\, X2) = det(Xo/ n — XiA±— X2A2) is the hyperbolic polynomial associated with A\,A<i. 
As was shown in Section HI the numerical measure ha is entirely determined by the eigenvalues \j(0), 
hence by the polynomial Pa- 

This in turn suggests a natural way to associate to any hyperbolic polynomial P of degree n in 
d+1 variables a probability measure hp ° n ^ rf > which coincides with the numerical measure \iA when 
d = 2 and P = Pa- Given a unit vector uj G S^ 1 C M d , let Ai(w), . . . , \ n (uj) G R be the roots of the 
polynomial equation P(A, u) = 0, and let -B^s) = B[Xi(uj), . . . , A n (cj)](s) be the normalized £?-spline 
with knots Ai(o;), . . . , A n (w). The "numerical measure" of P is then the unique probability measure 
Hp on W 1 whose Radon transform satisfies 

(TZ/ip)^, ds) = B aj (s)ds , 

where, by definition, (TZfip)^, I) = Hp({ x 6 M. d \ x ■ u G J}) for any interval I C R. 

In this generalized setting, the counterpart of normal matrices is the case where the polynomial P 
split into linear factors 

n 

P(X) = J](A - v k ■ (A l5 . . . , X d )) , with v k G M d . 

k=l 

When n > d + 1 and the vectors v k span the affine space M d , the density of hp with respect to the 
Lebesgue measure is the multivariate £>-spline in d variables, whose nodes are the v k, s. It is piecewise 
polynomial of degree n — d — 1. In the generic situation where any (d + l)-uplet of vectors v k is an 
affine basis, it is of class C n ~ d ~ 2 , see [6]. Again the density is log-concave in this case. 



39 



As in the two-dimensional case, the measure (ip can be expressed in terms of B u using the back- 
projection method. The example above suggests that, as the space dimension increases, the measure 
\ip becomes more singular. Then, the inversion formula has to be understood in the sense of distri- 
butions. In the three-dimensional case d = 3, we arrive at the simple expression 

HP = - ^2 A x J (x ■ uj) dcr(w) , (65) 

where da denotes the Euclidean measure on the unit sphere S 2 . As an example, in the particular 
situation where P = Xq — X 2 — X\ — X 2 , which corresponds to the differential operator d 2 — A of the 
wave equation, we obtain \ip = ^- da(x). This shows that, when d > 3, the support of fip does not 
need to be convex. Because the wave equation satisfies the Huyghens Principle, this example suggests 
that the link between polynomial regions of the density and lacunas of differential operators persists 
in higher dimensions. 
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